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Abstract 

The  aberration  of  plane  grating  monochromators  such  as  of  Ebert,  Czerny- 
Turner  and  Pfund  types  are  calculated  in  the  case  of  symmetric  arrangement. 
By  using  appropriate  focal  condition  and  grating  equation,  the  slit  form  for 
avoiding  the  wave-length  error  is  determined.  It  is  a  circular  arc  in  Ebert  and 
Pfund  types  and  is  an  elliptic  arc  in  non-concentric  Czerny-Turner  type.  The 
length  of  the  slit  restricted  by  astigmatism  and  the  optimum  size  of  the 
grating  derived  from  spherical  aberration  are  determined. 


1.  Introduction 

According  to  William  G.  Fastie,*’  both  the  entrance  and  exit  slits  in  Ebert  system 
are  to  be  curved  in  order  to  remove  the  effect  of  astigmatism  and  to  avoid  wavelength 
error  at  the  exit  slit  caused  by  the  curvature  of  spectral  lines  when  a  plane  grating 
is  used.  He  showed  that  the  resolving  power  obtained  with  any  short  portion  of  a 
long  slit  properly  curved  was  the  same  as  that  with  any  other  portion  of  the  slit.  In 
infrared  region,  in  which  radiation  energy  from  the  source  is  weak  and  detector  is 
not  sensitive,  slits  of  monochromator  must  be  long  in  order  to  let  much  light  energy 
enter.  Therefore,  for  an  infrared  spectrometer,  long  slits  are  very  desirable.  In  recent 
years  Ebert  type  spectrometers  with  long  circular-arc-shaped  slits  have  been  con¬ 
structed  by  many  workers*’  and  excellent  results  were  reported.  But  it  is  surprizing 
that  the  aberration  of  Ebert  type  monochromator  has  not  been  seriously  discussed, 
and  that  information  about  grating  equation,  focal  condition,  astigmatism  and  optimum 
grating  size  is  still  lacking. 

Since  Ebert  or  Pfund  type  monochromator  (Fig.  1)  incorporating  with  one  ccmcave 
mirror  as  collimator  is  a  special  case  of  Czemy-Tumer  type  which  has  two  concave 
mirrors,  we  shall  consider  here  only  the  two  mirror  system.  In  such  a  system,  be- 

1)  W.  G.  Fastie :  J.  Opt.  Soc.  Am.  42,  641  and  647  (1952). 

2)  For  example;  M. A.  Ford,  W. C.  Price  and  G. R.  Wilkinson:  J.  Sci.  Inst.  35,  55  (1957). 

D.  W.  Robinson :  J.  Opt.  Soc.  Am.  49.  966  (1959). 
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a)  Ebert  type  b)  Pfund  type  c)  Czerny-Turner  type 

Fig.  1.  Plane  grating  monochromator  using  spherical  mirrors  as  collimators. 


cause  of  numbers  of  adjustments  being  involved,  the  two  spherical  mirrors  are  apt  to 
be  so  placed  that  their  centers  of  curvature  are  not  coincident.  For  such  non- 
concentric  arrangement,  designers  of  monochromators  should  be  prepared  to  face  the 
problems  whether  the  curved  slit  is  applicable,  how  the  grating  equation  and  focal 
condition  are  affected  and  if  the  astigmatism  is  enhanced  or  the  higher  order  aberra¬ 
tion  can  vanish.  To  enlighten  them  on  these  subje:ts,  it  is  necessary  to  examine  the 
aberration  in  two  mirror  system  in  general.  For  this,  application  of  Beutler’s  theory” 
on  concave  grating  to  the  case  of  plane  grating  monochromator  will  serve. 

2.  General  Expression 

As  shown  in  Fig.  2,  the  origin  O  of  a  rectangular  coordinate  system  is  taken  at 
an 'arbitrary  point  on  a  holizontal  line  (x-axis)  which  passes  through  the  grating 
center,  and  z-axis  is  in  the  vertical  plane  and  parallel  to  the  grating  rul  lings.  Light 
from  a  point  (XiPiZi)  on  the  entrance  slit  Si  falls  on  the  point  XiXtVtZt)  of  the 
surface  of  mirror  Mi  and  reflected  to  the  point  G(f  v  1)  on  the  plane  grating  and  next 
diffracted  to  the  point  X'ix,'  p,'  z/)  on  the  mirror  M*  and  finally  reflected  to  the 
point  (x/  yi  z/)  on  the  image  plane.  Radii  of  curvatures  of  the  mirrors  Mi  and  A/* 
are  R  and  R'  respectively,  centers  of  the  curvature  being  at  Oi(xo  VoO)  and  O*(xo'  Vo'  O) 
respectively.  Grating  center  is  at  'COO),  its  width  w  and  the  angle  between  its 
surface  and  p-axis  0. 

Now  let  the  slit  Si,  the  mirror  Mi  and  the  grating  compose  Afi-system  and  the 
grating,  the  mirror  Afj  and  the  image  plane,  Af,-system.  If  we  take  the  point  Oi  as 
the  origin  for  Mi-system  by  parallel  transformation  of  the  coordinates,  we  have  for 
the  points  on  slit  Si,  mirror  Mi  and  grating  G 

for  slit,  Si=Xi— Xo,  Si=Vi— Vo,  Zi=Zi  \ 

for  mirror  x=x,— Xo,  V^V*— Vo,  z^z,  ^  (1) 

for  grating  Gi^f — Xo=Ao— ,  G,^v— Vo=«'^— Vo ,  /=/ 

3)  H.G.  Beutler:  J.  Opt.  Soc.  Am.  35,  113  (1945). 
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Fig.  2.  Coordinates  of  optical  elements  in  two  mirrors  system. 


where, 


ht=c—Xa,  0i=sin^,  ^tScosO 

Then  the  direction  consines  of  radius  R  of  mirror  Mi,  incident  light  path  p,  and 
reflected  light  path  pt  become 


X 

L=-^ 

Af=  ^ 

N= 

R 

,  X-Si 

R 

V-s, 

Vl  = 

R 

(2) 

Pi 

,  Gi-x 

Pi 

G,-v 

v,= 

Pi 

l—z 

Pi 

P. 

Pi  ; 
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and  therefore 


where, 


R^=x'+y'+z' 

Pi*  =  /?“+S*— 2(SiX+S,1/+2i2) 

P2*=i?*+{/*-2(G,x+G,y+/2) 

s2=S,»  +  S,*+£,* 
ff*=Gi*+G,*+/* 


(3) 

(4) 

(5) 

(6) 

(7) 


(8) 


s  and  g  are  respectively  the  distances  from  Oi  to  the  point  on  the  slit  Sj  and  to  the 
point  on  the  grating  G  as  shown  in  Fig.  2.  From  (1) 

(;*=(;o*+0fi+g»  (7') 

where, 

(;o*  =  Ao*+l/o* 

(/!=— 2(Mi+yoWM' 

go  is  the  distance  from  Oi  to  the  grating  center. 

Applying  the  reflection  law  on  the  surface  of  mirror  Mi,  we  have  for  the  direc¬ 
tion  cosines  of  reflected  ray 

vt=vi—2N{Li.\-\-Mt‘\+Nvi)  ' 

Using  (2)  (3)  (4)  and  (5>,  we  have  from  the  first  expression  of  (9) 

j^(PiG,-fftSi) 


(9) 


R*Pv—(h{p\'-s') 


Let  for  convenience 


Then, 


a=(Hlpi 

(10) 

ei=Pi*-^ 

(11) 

a=l—{aeilR*) 

(12) 

x={Gi+asi)a~^ 

(13) 

expression  of  (9) 

V={Gt+aSi)a-^ 

(14) 

z={l+azi)a-^ 

(15) 

and  hence  from  (3),  (6)  and  (7) 
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where, 

a*s*+2ak*+g*-R*a*=0 

(16) 

From  (1), 

5iGi + SjGj -|- 2i/ 

(17) 

where. 

**=*0*+*, 

(170 

ko*=hoSi  — VoSj 

(18) 

hi=i-Si^i+Sitl>t)w+Zil 

(19) 

Solving  (16)  for 

a  and  remembering  a  >  0,  we  obtain 

a:=-^  [—k*+Rsa(l  —  ma~*)^] 

(20) 

where. 

s*g*—k* 

R*s* 

(21) 

m  in  the  above  expression  is  obtained  from  (70  and  (170  or  (8)  and  (19)  as 

m  =  mo+mi+mt 

(22) 

mo= 

mi= 

m,3 


sW-V 

FPs'  ' 

jVi-V 

IPs* 

{(<4i*ro*+?>i*w,o+2^i^*5iS,)w*+  2(^>,Si— <4,s,)M;2i/+m,i/*} 


(23) 


»iio=Afl(s,*+2i*)+SiS,yo 
Wu = /loSiS, + (5i» + Zi*)yo 

ro*=s.i*+Zi* 

mti=Si*+St* 

With  the  value  of  x,  y  and  2,  we  obtain  from  (4),  (5),  (6),  (7)  and  (17) 

Pi*=/?*+s*-2(**+as*)a-* 

which  becomes  from  (20), 


(230 


Pi*=R'-+s*-2Rsa-ma-*)i 


(24) 
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Similarly  for  the  reflected  ray,  we  have  from  (5) 

Pi* = /?* -f- * — 2((;2 + aAf*)a  ■  * 

which  becomes  from  (20) 

DM 

(H'=R*+g^-2R'ma-^-2  (l-wfl-*)*  (25) 

Since  the  factor  a  in  the  expressions  (24)  and  (25)  is  an  unknown  quantity,  it  must 
be  found.  Using  (12),  (20),  (11)  and  (24),  we  can  obtain  a  relation  for  a  as 

0=1—  {— ^*+/?sfl(l— mo"*)*)  {/?— 25(1— mo"*)*} 

or 

(l+  *’-j-2»ia-‘+fl=  (/?*fl+2ife*)(l-»ia-*)i 

By  squaring  both  side  and  by  re-arranging 


/  /?*  \  A*  \ 

vKo)=  1-  ,  o‘+2  1-^  fl*-f 

\  s*  /  s*  / 


1-4  ^  (2** -f (7*1 


— 4mfl+4  =0 


Siuue  oasl,  approximate  value  of  a  is  obtained  from 


fl=l-y<l). 


'MV>^ 

da 


But  since 


V<l)=4(l-^-)(l-*’+m) 


do 


we  obtain 


0  =  1  + 


R* 


+  w*  + 


By  applying  (7')  and  (22)  to  the  above  expression,  a  can  be  expressed  as 

O  =  0o  +  0i  +  0i+03+04  + . 


where, 


flo=l+  3  ^0 »  />o=(Co+wio)Z?o 

Oi=  I  />i ,  />i=  Domi  -  (Co+nto) 


(26) 


(27) 
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_  2  ^ 

^1=  ■j'ft  » 

pi=Dom,—  {mi(2i+(7,(Co+mo)) 

2  ^ 

^*=  f 

P3=—  ^  (mxgt+mtQi) 

esileo 

III 

Q 

^  1 
p3=—  prm,gi 

C  =1  ^ 

4^2  . 

p 

III 

1 

“  1  *• 

These  Pi,  />,  and  />*  are  also  expressed  by  (8)  and  (23)  as  follows. 
2 

Pi—  jip  \{^iP\9'^^iPi\)n)—p\iZil] 

Pit—hf^Ctt—m^  —  Do 
Pii=yo(Co—mo) —  Do 
P»  =  -;rDo 

p,o=-(Co+mo)  +  4^^i«  +£>o  y 

— (Co+mo)+4-^j^"  +^0”?” 

2  s  s 

Pti=  j^g2  ihontii+yomio)+Do—^^- 

p20=  — (Co  +  /Wo)  +  /?o  ^1*' 


(28) 


(29) 


(30) 


+  i^pl'pn  +  +  ^1^^/>3«)W*2i/  +  (,<f>lpt1  +  ^ip3B)tOl*  +/>»*/*) 

/>3o=>«io+/io^o*  />3i  =  »»n+2AoSi5,+yo>'o* 

p3i=mio+homoo+2yoSiSt  p33=mn+yomto 

p3\=ho*-]-2hoSi  P3t=ko* — 22^0^] 
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/^3«  =  2(Si|/o  —  AoSj)  />37  =  WIio  +  /!oWIj1 

/>38  =  mii  +  yoW,i  Pxi»=Zika' 

Therefore, 

=  ^(1+^1+ ^1  +  ^3 +  64  + . ) 

6,  =  — 2flo~*fli.  6,=  — 2ao"‘o*+3ao"^ai2 
*3=  — 2ao"‘a3+6flo”^Oia*— 4ao"®ai* 

64=— 2ao"*a4+3(ao"*aj*+2ao"’'aia3'— 12ao"*ai*a»  +  5oo"*fli* 

a"‘=ao“‘(l+^i^  +  ^j^'i'^3^+^/  + . ) 

bi  =  2  2  — au"^ai‘) 

6s' =  2  {63-2a«-*ai(es-ai*)} 

64'=  2  {64— ao"*(a,*+2aias)+6ao'®a,-a,-3ao'‘a/) 

From  (22)  and  (32),  (24)  becomes 

Pi*=/?-H-s*— 2/?s  [^0— (wi+mo6i)— (W1  +  W161  +  W061)— (/w,6i+wi6,  +  Wo63) 

— (»i,6,+mi6s+»io64)- . (33') 

where 

^0=1— Woao"*  (34) 

Then  p,*  can  be  rewritten  by  expanding  the  last  term  of  (33')  in  series  as 


Pl*  =  PlO*  +  Pll'  +  Pls'  +  Pl3'  +  Pu'  + .  (35) 

p,o*=/^*+s*-2/^s^„i  (36) 

Pii'=/?sao"*^o"4(»ii  +  mo6i)  '  (37) 


p„'  =  /?sao"*^o‘^  m2+Wi6i  +  mo6j  +  -|-ao"*^o“'(»ii+Wo6i)*  (38) 

Pi3'H/?sao'*^o‘^  W26i+wi6,+  mo63  +  ^ao"*^o"*(»»i+wo6i) (»i,+»iA+Wu6,) 

+  ao-*^o"*(»ii + W06,)’  (39) 

p,4'  =  /?sao~‘‘^o"i(»i*6,+mi63+mo64)+  4j^fPii'*+2p„'Pi3')  (40) 

Using  (23),  (23'),  i29)  and  30),  we  have  p,/  and  p,,'  explicitly  in  terms  of  w  and  /, 
namely 


By  further  series  expansion,  incident  light  path  pi  is  finally  obtained  as 


Pl  =  PlO  +  Pll  +  Pll  +  PlJ  +  Pl4  + .  (41) 
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Pl*=  9  jpi/—  .  2  (2pii'p.3'-hPit'^)+~^  Pll'^Plt'—  c,  tPn*  (45) 

ipio  1  4pio  opio  O4pio 

Pio  is  the  path  length  for  the  incident  ray  to  reach  the  grating  center. 

In  an  exactly  similar  way,  reflected  light  path  ft  becomes 

ft* = + (7^ — 2/?*ao"M  Wo + (m  1  +  Wo^i') + (w* + m  A' + m  A') 

+{mtbi  +  w  A' + w  A')+ (wjfj,' + mibi  +  w  A') + . ) 

RkJ 

— 2  {^0— ao"^(wi + mof>i) —ao-\mt +mibi+ tn,M 

— ao"*(  w»6i + w  iftj  4-  m  A) — ao"*(  Wjftj + mifts + w  A) — . )  ^ 

By  series  expansion, 

ft* = fto*  4-Pi/4-ft»^+fts^4-ft«^+ .  (46) 

Rk  ^ 

fto*  =  /?*+go*-2/?*Woao-‘-2  $0^  (47) 

fti'=ai— 2f?*flo"Kwi+»iof>i'  +  ao‘*0o"^(wi+wofti)— 2  ^0^  (48) 

Dfc  2 

(h2'=gt—2R^ao~\m3+mibi'+mJ)i)+-  J* 

1  R 

m2+mibi+m^t+ -r +  —  ao"*^o'^^i(wi4-Wo6i)  (49) 
4  s 

p»3'  =  — 2/?*ao‘Hwi6i'+mA'+wA')+  Oo”*^o“^|w,6i+mA+Woft3 

+  i  ao”*^o”*(wi + Wo^i)  (w,  +  mif»i  +  Woftj)  +  Wi  +  »iof>i)® 

+  -y w,+wlA  +  wA+-jao‘*^o'*(wi  +  mA)*|  (50) 

ft4'  =  — 2/?*ao‘‘(w,6,'+wA'+Wo64')+  ao"*^o”^|^w,6,+mA+wA 

4-  ao‘*^o‘*  {2(mi  4-  Wof»i)  (nttbi  4-  wii6i  4-  w  A)  4-  ( Wj  4-  Wif>i  4-  w  A)*} 

3  5  1 

4-  -g  ao"‘^o'*(wi  4-  Wo6i)*  (mi 4-  mi&i  4-  m  A) 4-  flo**0o'*(mi 4-mo6i)M 

/?  I  1 

4-  —  ao~^$o~^  bi  I  mtbi  4- m ,fr,  4-  m A  4-  y  ao"*^o" ‘(mi 4-  mo&i) 


(m,  4-  m  i6i  4-  m  A)  4-  a„~*0o'Hm  i + ntobi)' 


(51) 
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From  (32),  (33)  and  further  from  (23)  and  (37),  we  obtain  for  p,,',  p,,'  and  pt/ 


p,/=g,(l-ao-‘)-y^oi(2*i+«.Pn')(/?- 

By  putting 

Pio  PjoS  \  S 
or 

/?-  ***  ao-’^o-*  =  (  ^  -Qo)ftosao-‘^o-* 
S  \  Pio  I 

Pti  becomes 


Pn'=—  *  Pii'-\-Bi-\-BiQo 

Pio 

fii=Pi(l-flo-‘)-2^“ 

Pio 

fii'  =  P»o«o*H2*i+aoPii') 


From  (8),  (19)  and  (37') 


(52 

(52') 


(53) 

(54) 


Bi-j-B/Qo—  — 2[  {?ii(Bio+fiio^Oo)+^(^n+^n'Qo)}M'+(Bij+fiii'Qo)^in 


B,o=hoa-a,-^}-  s.flo-* 

Pio 

fin=l/o(l-ao-‘)+  s,ao-‘ 

Pio 

D  —  /,  -1 

£fl2=  Oo 

Pio 


^io'  =  P*oPo~*|5i+  •'4ioj 

•fil/  =  P2oflo~*|  — S2+  j 

^12'  =  P20<3Io"'|  ~1+  •'^1*1 


(54') 


Using  (33),  (37)  and  (38),  we  have  for  P22' 


P22'=ff2-/?Wo-‘-  ao-‘^o-*) 

+  Pii'*+  Pii'+/?*»»oao**ai* 
Further  from  the  relations  (23)  and  (52') 


P22'= ^^Pi2^  +  B2  +  Bi'Oo 

Pio 

B,=P,(l-ao-‘)+  ^  Pii'j  +  |-/?*ao- Wi*  I 

=  Pis'p2oao~*  1 


(55) 

(56) 


From  (8).  (19),  (29),  (37')  and  (33'), 
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Bi + Bi  Qo  —  {4>i\Bit+ Bio'Qo) + ^zKBii + B21  Qt)+^i^^2Bti + w* 

+  {^i( — 2fi2j+fijj^Qo)+9^(2fij4+3»/Qo)}  wzj  -\-{Bii+ Bti'Qo)l^ 

B2.=l-flo-‘+  >ln) 

B22=-  ^10) (s2-  ^^5^")  +  "^  moMn 

. 

D  _  «o"‘  / 1  Oo  .  \(  Oo  -  \  16ao'*  „  .  ^ 

B2,  =  l-a,-^+  (1-  moA2*2i* 

Bio'^PioOo  *-420  »  B2i'  =  P2oOo  *^21  »  ^22  =^20^0  *-<42* 

•fi23^  =  PtO®o”*-'‘l2S  •  ft/  =  Pj0®0  *''4i4  »  •Sis  =  PlO^o  '.Ali 

The  second  term  of  #>«*'  (50)  becomes  (Ar*/s*)  by  the  use  of  (39).  Therefore  from 
(28),  (32)  and  (33),  pzs'  is  finally  expressed  as  follows. 

P2i=~  Pi/  +  Ss  +  P2oPiS^Qo  (57) 

Pie 

Be=  Pn'Pi2  +  ll  P.2'  +  |-/?^flo-»m.^.* 

+  |-/?*Woao‘*/>i/>2--p-/?*»»oflo'*/>i*  (58) 

or  from  (19),  (23').  (29),  (30),  (37')  arid  (38') 

Ss = ((4i*Sso +^i*<^S3i + ^\^Bn + 

+(^i*Sj4+^i*S»s+<^i^iSM)M;*Zi/+(^iSs7+^ftg)M;/*+Ss»2i/* 

€1  1  32 

Sjo=—  9/f«s*  ®®~*^*®^*®* 

.  16  ^  128  , 

+  WoOo  */>io/>2o--^»»o«« 

Sji= - (./1iO'4i1"1‘-'4ii./42o) - pp(Sl-A22~S2.42o)+ . 

Bi2=  ~  (Si.421~S2-'4i2) 
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32ao”*  16 

“  9/^451  (^•®/^l»*  +  2»li|^10^ll)+ moflo  *(/>l0^l+/>ll/>l2) 

128 


“  "jljp  ^11^*1  +  ■^Sj/ln  +  " 

Bu  =  — (i4  ti^io + A 10'^ m) + • 


Bii=—  (/I  i2>lji+i4 11^424) +• 


Bit=—  >1 12^28  +  -^  >^26  —  ht'pxi'Zx* 

~  ^o/>12/>2S+4  ^  ntop\2*Zl* 


Similarly, 

P2t'— - Pl/  +  ^4  +  P2oPl4^0o 

PlO 

®4=  ■^^(Pi2^*+2pi/pi80+  Pis^  +  'g  R*ao~^mipx*+  g’/?’oo”*^i/>i/>2 

-  y/?»moOo'*(/>2*+2/>i  R*m^o~*P\'p2 

By  further  series  expansion,  P2  is  finally  presented  as  follows 


P2  —  P20  + P21  +  P22  +  P23  +  P24  +  • 


(H\  =  "IS - P21 

2P20 


Therefore,  optical  path  traversed  through  A/|-system  is 


1 
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Pi+fti— (P10+P20I+ ■«  r  „  Pi/+  „  P2/I 
^  L  Pio  P20  A 

,  1  r  1  ^  '  1  1  '.1  • 

2  L  Pio  P20  4pio*  4p2o*  J 

.  1  r  1  , .  1  ,  1  ,/  1  , .  1  ,  1  1  ,.\ 

+  ^  Pis  +  .  P2S  —  n  ,  P2I  I  Pit  +  P22  —  j  Pll  *—  .  ,  Ptl  * 

2  L  Pio  P20  2p2o*  \  Pio  P20  4pio®  4p2o*  / 

—  P12  I  *  Pll  “  „  2  P21  +  a»  s  1  „  •  P"  ~  -  *  P21  I  I 

2pio  \  Pio*  P20*  /  oPio*  \  Pio"  P20*  /  J 

1  r  1  /  .  1  .  1 

2  L  Pio  P20  2p2o 


J_„  ^ .  _L„  ^ 

Pis  -f  0  P2S 
PlO  P2" 


7P*i 


il/i  s  Pn'*—  A  3  P2/*]—  o  (pi2*“  A  ,  Pi/*) 
\  Pio  P20  4pio*  4p2o*  /  2pio  \  4pio*  / 


/  *  « ' 

1 

n  '  1  ^  n  '  n  '' 

1  p,.*  ■ 

P..-  ) 

2pio 

Pis  +  ,  Pll  Pl2 

2pio® 

P..'  ■ 

■  P»'  ) 

+  si  pI.‘  si'  '^■'‘)]+ . 

Inserting  (53),  (55),  (57)  and  (59)  into  the  above  expression, 

Pl  +  P2  =  (PlO  +  P2o]  +  '^j^"^  (fil  +  ^/Qo)J 

+  i  [1  <«■  +  *'«?.)-  +fi.'<?.r] 

(pi. ‘'"'■■I'®')  hall'"'" 


1 

L  .  M 

-L 

PlO 

+ 

PlO  P20  / 

“  P..' 

. 1 

■  2lr(‘’'’'+  w  . ) 


1 

^  .  /!_  P 

1  1 

1 

8 

J  Pll  Pl2  - 
PlO 

P20*  ' 

1  PlO 

1  , 

1 

8p20* 

,  1 
P20* 

Pll'*-I- 


. ro 


(66) 
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Pio  is  the  path  length  of  reflected  ray  to  reach  the  grating  center.  As  will  be 
mentioned  lafer,  Qo=0  (52)  is  approximately  the  condition  for  image  formation,  there¬ 
fore  Qo  is  made  zero  in  third  and  fourth  terms.  Furthermore,  from  (37'),  (38')  (54'), 
(56')  and  (58'),  (pi+Pi)  is  explicitly  expressed  in  terms  of  tv  and  /  as  follows. 

Pi  +  Pi  =  [PlO  +  Piol  +  ~2  +  9^1 1)**' +^12Zl/] 

■F  ~2  +  +  ^26/*l 

+  [(<4iVm  +  9^1*?^S1  +^1<^V'S2  +  ^*V’Ss)M'* 

+  (^1  Vm  +  ^  +  (<4|^S7  +  ^^s»)iol*  +  v'’s*^lf  *1 

+  \[Mu’*+l'y+ . ]  + .  (66') 

For  A/2-system,  an  analogous  expression  for  {p/ +Pt')  may  be  derived  by  the  following 
transformation  of  coordinates. 


So  — = 

0 

1 

— V , 

/?-/?' 

Pu'-*Pu''.  P12'  — Pll", 

Pis'  —  Pu"  , 

Pu'  -*  Pu' 

F.-F,,  Fi'-F,', 

F2-F0, 

Therefore,  the  total  light  path  length  or  optical  path  function  L  is 

r  ,  ,  ni. 

^=Pl-fP2  +  P2  +Pl 


(67) 


where  n  is  the  order  number  of  spectrum,  d  the  grating  constant. 

If  L  is  constant  independently  of  w  and  /,  the  Fermat’s  principle  will 
That  is. 


3/ 


=0 


be  satisfied. 

(68) 


This  means  that  the  optical  system  is  free  from  aberration.  However,  it  is  generally 
difficult  to  satisfy  the  above  two  conditions  at  the  same  time,  therefore,  the  system 
is  not  wholly  devoid  of  aberration.  The  residual  aberration,  however,  can  be  made 
considerably  small  by  proper  alignment  of  optical  parts.  To  deal  with  this  aberra¬ 
tions,  it  is  convenient  to  devide  it  into  various  known  aberrations. 

Therefore,  following  Beutler’s  method,  let 


L=  F\-\-Fi +(^i+'F2')+(Fj-f Fj')-f (^4-)- F4')-1- . 

d 


(69) 


'  Fl  =  Pl0-fP20+'^(^l<('l0  +  ?^l^^u)W 


(70) 
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=  y  (9^1  V*0  +  0*^21 + 


+  (9^iV»o+9ii*<M'’3i  +  <ii9i2Vs2+9i2Vss)w* 


/^3=  2  +  ~2 


1 


Fi=  2<^io(u>'+l*)' 


(71) 

(72) 

(73) 


where.  0io,  v^u,  .  are  expressed  as  follows  by  referring  to  (66),  (66')  and  from 

(37'),  (38'),  (54'),  (56'j  and  (58') 

2 


010= —  ^  (-^lO  +  filO^Qo) 

P20 


011  = -  (fiii  +  5i/Oo) 

P20 

012  — - (^12  +  fil2^0o) 

P20 

020=  (^20  +  -fi2n’'Qo) —  jn,  ,  •^lO*" 

P20  Pl0*/?*5*  P20 


021=  (^21  +  ^2/Qo)~  SPlc* 

P20  Pio*/?*5* 


""  pi  ~  PioW  ‘  • 

P20  (^20  +  fi2/Qo)-  p,„S/?252  >1i2*2i* 


(filo'f  ^lO^Qo)  +  p  >ll 
P20  Pio«S 


(74) 

(75) 

(76) 
'(77) 


1 

1 

P20 

P20 

(Bij-t-Bii'(?o)+ 


PioFs 


Au\  (78) 


1 

P20 


(^12  +  fil2^(?o)+  p 

P20  Pio/rS 


«1* 


,  _  1  P  1  ,  A  to  B\o  \  ,  1  /  4  ^10*  \ 

0SO=  ^30 - 02ol  r>  — ^  '  1  "!■  I  ■^20 —  ,x>*  1  ) 

P20  P20  \  PioFs  Pzo  I  PlO  \  PlO*^S*  / 


■^10  /  1  \  _  •filO 

Piof?S  \  PlO  P20  /  P20* 


031=  531  +  ’ 

P20 


(79) 

(80) 

(81) 

(82) 


,  _  1  P  1  r  (  Aio  Bto  \  ,  I  All 

9^32=  B32—  0’21  ‘  D  “  +022  ■  p 

P20  P20  \  Pl0«5  P20  /  \  PlO*? 


fill 
Piofis  P20 


*4io  / 

~+  — 1- 

filO  1 

P,0*fi*S*)l 

Piofis  \ 

PlO  P20  / 

P20*  1 

4-(a 

_J®‘L 

/ 1  PioBs 

\  PlO  P20  / 

P20* 

(83) 
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*l>t\  '  ■^u 
P20  PioRs 


Bu  \ 

PiO  I 


+ 


^f’s»=  _  OS9 - D-  ~ 

P20  P20  Pio«S  P20  I 


(84) 

(85) 


Lastly,  it  is  clear  from  (66)  that  the  term  which  includes  =(«;*+/*)*  in  ipw  is  only 
(1/Pjo)  Bi,  therefore  the  term  which  includes  only  in  (60)  is 

Pi2'*+|^/?*WoOo->.* 

Further,  from  (38)  and  (28)  the  above  becomes 


{/?s/7o'*^o”^(»»2  +  WoA2))*+^  /?*moffo'*|7)om2-  ^(Co+»«o)| 

Since  the  terms  which  have  in  mo  (23)  and  60  (32)  are  ga//?*  and  2ao~^ai^  =  ■|•ao“'/>2 
final  expression  for  ^40  is 


3 is'  %  ’^)1 

Analogous  expressions  for  Fi',Ft'  etc.  can  be  obtained. 


(86) 


3.  Symmetric  properties  of  monochromator 


Let  us  consiner  the  case  in  which  «=0  and  9=0,  that  is,  ?>i=sin^=0  and  ^= 
cos^=l.  This  case  corresponds  to  the  zero  order  spectrum  or  plane  mirror-like 
reflection.  In  such  a  case,  most  terms  of  the  optical  path  functicMi  L  vanish.  For 
such  setting,  consider  the  condition  for  the  first  order  aberration  to  become  negligible. 
First,  the  equation  (70)  becomes 


and 


Ft=Pio+Pio+Y*^^'^ 

F/=p,o+  P20'  -b  y  V’l 


When  Fermat's  principle  is  applied  to  (Fi+Fi), 

=0 

Therefore,  from  (75)  and  (54) 

(l-flo-‘)+— (l-ao'->)  +  -^flo' 
P20  Pit  Pto  Pit 


-ao-(s2-  *~  An]Qo-a„'-Hs,'-^At/  1<?,'=0 
i\s  \  tc  9 


(87) 
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The  conditions  that  satisfy  the  above  equation  are 

yo=-KVo,  Si'  =  KSi,  s/=-«S2,  ] 

>  (88) 

Zi'=—KZi,  R'  =  kR,  ho=Kho' ,  ) 

and  therefore,  from  s(6),  flfo(8),  ^0(18),  »io(23),  flo(28),  ^9(34),  Pio(36),  P2o(36),  mii(23'). 
AuiST)  and  Po(52) 

s'  =  KS,  Pio'  =  *Piot  P2o'  =  *P2o,  ao'=ao,  i4u'  =  — *Mii,  Qo'  =  —k'^Qo  (88') 
where,  *>0. 

Since  the  case  of  Ar^l  has  no  practical  meaning,  only  the  case  of  *=1  will  be 
considered.  Then,  (88)  is  the  conditions  that  Mi  and  M2-systems  are  symmetric. 
One  of  the  optical  arrangements  satisfying  the  above  conditions  is  shown  in  Fig.  3. 


Fig.  3.  Schematic  diagram  of  symmetric  optical  arrangement 
of  two  mirrors  system. 

In  such  a  case,  even  though  or  0=^0,  the  conditions  of  (88)  remain  valid,  and 
^22+(/>2i,  V’si+V’s/.  <Pia+</>as'>  •pi2Zi+(f>i2Zi'  and  (/^^tZi+tpstZi'  all  vanish  because  of  their 
asymmetry.  Therefore,  it  is  sufficient  to  consider  only  the  following  terms  instead 
of  (70)~(73). 
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•f'i+^/=2(pio+P2o)+^j<^ioW/  (70') 

Ff\-Fz  =  (^4i*^2o + *^so + ¥^1^  (71 ' ) 

F,+F,'=<^.5P  (72') 

F^+Fi'=iPUu}^+P)'  (73') 


Although  aberration  expressions  for  Czerry -Turner  type  monochromator  in  general 
case  are  too  complex  to  examine,  if  the  system  becomes  symmetric,  those  expressions 
are  simplified  and  we  can  easly  derive  various  conditions  for  image  formation. 


4.  Focal  condition 


Since  the  grating  equation  in  the  case  of  long  slits  can  not  be  derived  unless  the 
focal  condition  is  given,  we  shall  consider  it  for  the  image  formation  in  the  first  place. 
This  should  be  determined  from  o(Ft+Ft')ldw.  Therefore,  from  (71') 

(¥>i  V:o + <4= VsOm' +((i.  Vm + = 0  (89) 

The  solution  that  always  satisfies  the  above  equation  independently  of  w  and  0 
should  be 

=  —  —  ^  (90) 

Inserting  Ai(37'),  ^2(88'),  jBi(54'),  ^2(56')  and  Bs(58')  into  (77),  (78),  (81)  and  (83), 
furthermore,  using  mi(23'),  p\{2S)  and  />2(30),  we  see  that,  if  Qo+^=0,  (90)  is  satisfied, 
and  since  J  is  negligibly  small,  the  focal  condition  is  approximately  given  by 

Q.=  ^ - )=0  (52') 

Pio  PioS  \  S  j 


This  equation  is  very  simple  in  form,  but  it  is  not  so  easy  to  obtain  s  or  Si  from  it. 
However,  approximate  solution  is  obtained  as  follows. 

Since 


“+  — 

PlO  P2oS 


Rao0o^ 


multiplication  of  (52')  and  the  above  expression  gives 


¥>(S)=P20*S»-p,o*(/?flo^O-  -  )  =0 


Now,  from  (28),  (34)  and  (23) 


1  2  ‘  1  1 

= 1  -  2  ^0+ y  y  -g  mo*+ 


:1- 


2R*S* 


+  ■^1  Co+ 


. 


(91) 


(92) 
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The  values  given  in  (28)  and  (34)  make  pio*  (36)  as 


D  =  1  — 
£?o-l  ^ 


Similarly, 


Vn  /  4  1  4  \ 

Pl0*  =  (^  — «)*+  . I 

P2,.=(/?-  j )(i-f/>®+  9-/»o’+ . ) 


4-2^® 
+  35* 


(l-|/>®+ 


Therefore,  the  equation  (91)  is  expressed  as  follows. 


,<»)=- |s«-2,)+^l  I  («- 

-  -'r)  (‘'"+  ) 


2/?s*  3  ^^®  ^®+  Rts* 


Since  s»/?/2,  approximate  values  of  s  and  Si  «u:e 


ro*=S2H2i* 

> 

The  following  calculations  are  necessary  to  obtain  the  value  of  <p(RI2) 

(  V^o)»  >•  */»  ~  ^0*  ^9*  +  To 

where, 


Y,=v>  <•.')* 
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On  the  assumption  that  Ao,  the  x-coordinate  of  the  grating  center,  is  near  /?/2  as  in 
customary  arrangement,  1— 4Z)o/3  becomes  negligibly  small. 

Hence 

v’(^) =-2(AoVoH  yo)|  (  a:.+ •Jri/os.jV  |-z>o  v,s.  j  [ 

. 

=  —2(ho^ro^+Yo)^Ki-{--^yoSy^  |l+  ^  Z)o 


Further,  the  following  calculations  are  needed  for  obtaining  d<p{RI2)lds. 

i-i 


from  which 


(  Vo^)f-JI/J  — Po 


2hoS 


-hH 


-i 


+PVa 


=  y, 


whence  we  have 

MP/2) 


ds 


=  2/?» 


Ki+  1/0S2J 


*  ,  Wr,^+Y,) 

I  r> 


+|o.(i-|o.)]-2r,[(/f,+^»,s,)’ 

+  "^^o|^i+  yoS2j^+ . 

=  2/?>(a',+  ^  yoS2)*[l+  ^  (hoW+Yo) 

. ] 


which  leads  to 


R  ,  Ao"rb*+Fo 
2  /?>  ' 


(99) 


1  +  "y  VoSi 

F, 


)“l[i-Jr(*oV+F.) 


/?>  3  pzVoSi  . ] 


(100) 


( (101) 


sic={sc^-rt^'>- 
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where,  suffix  C  signifies  the  non-concentric  Czerny-Tumer  type  (i/e^O). 

As  for  Ebert  or  Pfund  type,  since  yo=0  and  therefore  since  Yt=Y\=(i,  we  have, 
with  suffix  E  to  signify  this  type, 


Sk  — 


-^4. 

hoW 

2  ^ 

/?* 

1+  “ 


R* 


hoW+- 


5iB={Sa^-ro^)^ 


(102) 


(101)  and  (102)  lead  to  the  conclusion  that 

(1)  In  the  case  of  non-concentric  Czemy-Tumer  type  sc  and  Sic  change  as  the 
coordinates  (s>2i)  of  object  point  change,  which  means  that,  there  is  no  fixed 
focal  plane  on  which  x= const. 

(2)  In  the  case  of  Ebert  or  Pfund  type,  if  ro= const.,  it  becomes  s«= const,  and 
Si <»= const.  In  other  words,  when  object  points  are  on  the  circle  of  radius  r^=yi), 
unage  point  are  also  on  the  same  circle,  symmetrically  opposite  to  the  former. 

Since  the  grating  center  ho  and  the  midpoint  of  slits  can  be  at  any  place 
unless  the  portion  on  the  collimator  surface  to  be  used  is  de^gnated,  we  can  easily 
determine  from  (101)  Sc®  or  Sic®  for  the  midpoint  S20  of  slit  and  from  (102)  s«  or  Si^ 
for  all  image  points  as  shown  in  Fig.  4.  When  the  portion  on  the  collimator  surface 
to  be  used  is  designated,  S20  is  uniquely  fixed  if  ho  is  given. 


Fig.  4.  Schematic  diagram  showing  the  relation  of  entrance  slit  and  image 
plane  in  symmetric  optical  system. 


But  when  in  Ebert  type  ro^  const. -the  case  of  straight  slit  for  example —focal 
points  will  not  be  on  the  plane  on  which  x= const.  In  a  symmetric  system,  the 
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image  of  straight  slit  is  apt  to  be  considered  also  straight  in  form.  This  may  happen 
in  the  case  of  extremely  short  slit  for  example,  but  not  necessarily  so  with  a  long 
slit  because  of  the  off-focus.  In  a  monochromator,  the  form  of  slits  must  be  de¬ 
termined  to  minimize  the  wavelength  error  arising  from  the  turning  of  grating.  For 
this  determination,  the  grating  equation  should  be  discussed. 


5.  Grating  equation  and  image  form 

Using  (70'),  (74),  (54')  and  focal  condition  (52'),  we  obtain  for  Fermat’s  principle 

Since  the  grating  in  Fig.  2  is  used  in  negative  order,  let  w  be  — Then,  from  (36) 
and  (47),  we  have  for  the  grating  equation 

«i=2dsin».//  (103) 


where 


H  =  a„->  1  + 

s  s* 


h^ao-'-l) 
-  R-Oo 


go-  ]“* 

(R-gor 


(104) 


a)  Case  of  Ebert  or  Pfund  type. 

We  have  seen  that,  in  the  case  of  Ebert  or  Pfund  type,  if  Si  and  Tq  are  given 
and  therefore  if  5  is  given,  the  focal  condition  is  satisfied  for  all  the  image  points. 
Since  ko*=hoSi,  and  both  Oo  and  Po  are  functions  of  only  s,  H  is  a  function  of  s  and 
Si.  Therefore,  H  becomes  fixed.  That  is,  the  incident  ray  of  wavelength  i  frqm 
any  point  on  a  circular  arc  of  radius  r^=yi)  with  its  center  yo=0  is  imaged  without 
wavelength  error  on  the  symmetrically  opposite  arc  of  the  same  circle,  and  the  wave¬ 
length  on  the  plane  of  image  circle  is  changed  only  by  the  turning  of  grating.  This 
image  circle  is  none  other  than  the  well-known  “  Ebert  circle  ”  considered  by  Fastie. 
Therefore,  Ebert  type  has  a  strong  point;  there  is  no  wavelength  error  even  if  long 
slits  are  used. 

The  relation  between  object  and  image  in  Ebert  type  can  be  explained  by  Fig.  5. 
Let  the  grating  G  be  a  circle  in  form.  In  the  figure,  the  outer  circle  Ci  is  the  colli¬ 
mating  mirror,  E  is  the  slit  circle  or  Ebert  circle  and  Cj  is  the  intersecting  circle  of 
the  collimating  mirror  with  the  sphere  whose  radius  is  pzo  and  center  is  coincident 
with  the  grating  center.  Since  the  circle  C*  is  uniquely  determined  when  Ebert 
circle  is  given,  let  it  be  named  “  Sub-Ebert  circle.”  M  and  M  are  portions  on  the 
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collimating  mirror  surface,  the  sizes  of  which  correspond  to  that  of  the  grating  and 
their  centers  are  on  the  sub-Ebert  circle  and,  further,  they  are  symmetric  with  re¬ 
spect  to  the  grating  center. 


Fig.  5.  Schematic  diagram  of  image  formation  of  Ebert  type  monochromator. 

C„  collimator;  E,  Ebert  circle;  Cj,  sub-Ebert  circle.  / 

0 

Now,  the  light  from  point  A  on  the  circle  E  is  reflected  to  the  grating  G  by 
the  mirror  portion  Af,  and  the  diffracted  ray  Is  imaged  on  the  point  A'  on  the  circle 
E  by  the  mirror  portion  M'.  If  the  object  point  A  moves  along  Ebert  circle  E,  Af 
and  Af'  move  always  along  the  sub-Ebert  circle  C2  and  the  image  point  A'  moves 
always  along  the  circle  E.  Since  A'  and  Af'  are  always  symmetric  to  A  and  Af  with 
respect  to  grating  center,  aberration  by  Af,  if  any,  is  compensated.  The  figure  shows 
also  that  the  maximum  movable  range  of  Af  and  Af '  along  the  “  sub-Ebert  circle  ”  is 
the  range  in  which  the  collimator  surface  is  eff^ectively  used. 

For  Pfund  tjrpe,  because  of  yi  being  very  small,  the  same  as  mentioned  above  is 
applicable  provided  that  the  slit  is  not  long. 

In  the  case  of  straight  entrance  slit,  values  of  s,Si,  and  ko  vary  along  the  slit. 

Hence  //  has  various  values  despite  the  angle  d  of  the  grating  being  fixed,  which 
means  that  the  wavelength  of  diffracted  ray  on  the  point  symmetric  to  the  object 
point  differs  from  that  of  the  image  center.  Thus,  in  order  to  minimize  the  wave¬ 
length  error,  form  of  the  exit  slit  is  bound  to  be  made  to  differ  from  that  of  the 
entrance  slit.  Such  a  case  is  not  a  symmetric  system  in  result.  In  such  an  asym¬ 
metric  system,  (60)  that  includes  cos  0  should  be  looked  into.  Since  the  focal  condition 
for  a  point  on  the  image  is  different  from  that  for  any  other  point,  a  definite  image 
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plane  can  not  be  fixed,  and  the  inclusion  of  cot  0  in  H  can  only  be  understood  if  it 
is  on  the  assumption  that  every  image  point  lies  on  the  focal  plane  for  the  center  of 
the  image.  Therefore,  in  the  case  of  a  straight  entrance  slit,  the  curve  for  the  same 
wavelength,  the  “  isochromat ",  will  vary  in  form  as  6  is  varied  which  fact  will 
deprive  the  monochromator  of  its  proper  function.  So  we  shall  consider  only  for  the 
symmetric  system  the  conditions  under  which  the  entrance  and  exit  slits  of  the  same 
form  can  be  used  and  H  becomes  a  constant  independent  of  the  form  of  the  slits. 

b)  Case  of  the  non -concentric  Czemy-Turner  type 

If  i/o  is  negligibly  small,  the  case  may  be  considered  to  be  that  of  Ebert  type. 
But  if  Vo  is  too  large  to  be  neglected,  all  the  conditions  for  image  formation  are 
affected  by  yo-  For  the  exit  slit  to  give  pure  monochromatic  ray,  it  should  be  placed 
along  the  isochromat.  It  is  therefore  convenient  to  have  the  form  of  the  isochromat 
to  be  independent  of  the  wavelength,  for  the  slit  form  remains  set.  For  this,  H 
should  be  independent  of  the  coordinates  of  the  slits  and  remain  as  a  constant  which 
is  Ho  for  the  midpoint  (sio,  Sio,0)  of  the  slit.  As  mentioned  in  the  previous  section 
concerning  the  focal  condition,  Sio  is  directly  obtained  from  (101)  when  sjo  is  given. 
Thus,  along  the  isochromat 

Ho=H 

From  (101)  and  (104),  H  expressed  by  any  point  (si  s-.  Zi)  of  the  isochromat  is  a  func¬ 
tion  of  SjZi.  That  is,  it  must  be  Ho=H{s2Zt)  Accordingly,  put 

9{s-,2i)=0=Ho-H{S2Zi)  (105) 

From  (104)  it  is  clear  that  predominant  term  in  the  grating  equation  is  sjs  and  the 
effect  of  Vo  to  the  value  of  H  is  very  small.  Therefore,  the  difference  between  Zi  of 
non -concentric  Czerny-Turner  type  and  Zih  of  Ebert  type  at  the  same  Sj  will  be  very 
small,  provided  that  ho  and  yxo  are  the  same  in  both  types.  Thus  Zi  becomes  ap¬ 
proximately 

^  .  f5#(S22i/») 

2i  =  Ziii  —  0{S<Zib)i  3_  " 

/  OZi 

=ZiB+  {//o— //(SoZib)}/  {3Ms2Zi/»)/3zi}  (106) 

From  (101),  (23),  (28),  (34)  and  (98) 


(107) 
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Inserting  a  given  5«  and  also  Ziw  at  s*  into  (101),  (104),  (107)  and  (103).  we  obtain  Zi 
from  (106)  and  when  those  s-i  and  Zi  are  used,  s  or  on  isochromat  is  determined 
from  (101).  For  example,  results  of  numerical  evaluation  are  shown  in  Figs.  4  and  6. 
From  these  figures,  we  see  that,  if  yo<0,  the  isochromat  resembles  an  oblate  ellipse 
and  always  Si  ^  Sio,  and  if  i^o  >  0,  it  resembles  a  prolate  ellipse  and  always  Si  ^  Sio. 

If,  in  Fig.  5,  we  replace  the  collimator  mirror  by  two  collimator  mirrors,  Ebert 
circle  by  entrance  and  exit  semi-elliptic  arcs  and  sub-Ebert  circle  by  entrance  and 
exit  sub-semi-elliptic  arcs,  it  will  very  well  demonstrate  how  the  image  is  formed  in 
a  non-concentric  system.  But  this  can  not  be  tried  in  the  case  of  long  slits,  for  it  is 
difficult  to  construct  slits  of  such  shape.  Therefore,  if  for  some  reason,  non-concentric 
Czerny-Tumer  type  has  to  be  used,  Vo  must  be  made  as  small  as  possible,  or  else 
the  slits  must  be  short  as  possible  as  in  the  case  of  Littrow  type. 


6.  Astigmatism 

Astigmatism  arises  parallel  to  the  grooves  of  grating  and  is  small  when  the  slit 
is  short.  Therefore,  in  the  case  of  short  slit,  even  if  it  is  straight,  the  spectral  re¬ 
solution  is  not  very  much  affected  by  astigmatism.  However,  as  the  slit  becomes 
long,  astigmatism  is  enhanced ;  therefore,  even  if  the  slit  is  of  circular  or  elliptic  arc, 
it  must  not  be  too  long.  For  accurate  determination  of  astigmatism,  ray-tracing  is 
needed,  nevertheless  we  shall  calculate  it  roughly  as  follows.  Assuming  that  the 
focal  plane  is  approximately  perpendicular  to  pio',  we  obtain  approximately  for  the 
aberration  Jz  from  the  image  point  z'(=— z)  by  the  central  ray  as 


Jz= pio'  (Fj-f-F»') 

Since  our  optical  system  is  symmetric,  from  (72') 

JZ  =  //2|/| 

where, 

H2  =  2picfp2i 

From  si>M(80),  B2t(56'),  Ai2(37')  and  Bi2(29),  H2  is  given  by 


(108') 


(109) 

(109') 


“-moW] 


PioRs 


9p2oR*s* 


(111) 


where, 


A 12= Oft  *^o  ^  ho^  1 —  Q  ntt)Do 
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Only  the  component  Jz'  which  is  perpendicular  to  isochromat  affects  the  spectral 
resolution,  and  this  component  is 


where, 


To'  is  the  distance  between  the  point  (sa  Zi)  and  the  point  at  which  the  normal  to  the 
isochromat  at  the  point  (ssZi)  intersects  y-axis  (Fig.  6).  The  isochromats  in  the 
range  of  Zi  given  in  the  figure  are  all  approximately  expressed  by  circular  arcs,  and 
from  the  figure  we  obtain  ro'=13.5cm  for  2/0 =5 cm,  and  ro'=7.5cm  for  yo=— 5cm. 
The  result  of  numerical  calculation  of  H3  for  the  setting  of  Fig.  4  are  shown  in 
Fig.  7. 


Aberration  coefficient  H,(cra) 


Fig.  7.  Effect  of  Astigmatism  on  spectral  resolution. 


From  the  figure,  we  see  that  the  effect  of  astigmatism  on  the  spectral  resolution 
becomes  gradually  large  as  the  slit  length  increases  even  if  circular-arc-shaped  slit 
(Ebert  type)  is  used.  Thus,  even  in  Ebert  type,  slit  length  is  restricted  by  astigma¬ 
tism  which  can  not  be  corrected.  Furthermore,  the  figure  shows  that  Jz'  increase  in 
the  order  of  the  cases  yo=5,  i/o=0,  2/0=— 5.  Since  Fig.  7  is  the  case  of  yio=const., 
as  seen  from  the  values  of  //  (Fig.  4),  the  portions  of  mirror  surface  to  be  used 
diverge  from  x-axis  in  the  order  of  the  cases  yo=5,  0,  —5.  That  is,  the  degree  of 
off-axis  of  the  effective  mirror  surface  with  respect  to  the  grating  increases  gradually. 
Therefore,  such  a  large  off -axial  arrangement  as  1/0=— 5  is  not  commendable. 
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'7.  Spherical  abwration  and  optimum  grating  size 

The  quantity  representing  the  spherical  aberration  is  (F4+F4').  Therefore,  in 
the  symmetric  system,  from  (73') 

F4  +  F4'=04o(wH/*)* 

^40  is  given  by  the  equation  (86)  and  is  independent  of  d.  Following  Beutler,  let  us 
consider  a  circular  grating  of  radius  .  The  condition  for  image  formation 

is 

Therefore, 


Since  the  effect  of  po  is  small,  only  the  case  of  Ebert  type  will  be  calculated.  In  the 
case  of  the  setting  shown  in  Fig.  4,  optimum  diameter  2a  of  the  grating  is 

5.392  {i(/i))i 

2a  is  shown  in  Fig.  8.  If  ho~RI2,  <p^a  becomes  approximately  (1/2F*). 

Therefore, 


that  is,  the  optimum  grating  size  is  approximately  proportional  to  R^  provided  that 
the  grating  is  ideal,  and  if  not,  the  size  is  reduced. 


Wavelength  in  /<. 

Fig.  8.  Optimum  diameter  of  grating  in  Ebert  type  monochromator 
in  the  case  of  /?=100cm,  A„=50cm,  yo=lOcm, 
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8.  Conclusion 

From  the  foregoing,  the  following  conclusion  was  drawn  concerning  the  plane 
grating  monochromator  regarded  as  of  perfect  symmetric  system. 

a)  With  Ebert  and  concentric  Czerny -Turner  type  monochromators  provided 
with  circular-arc -shaped  slits,  the  wavelength  error  can  be  removed  even  if  the  slit 
are  fairly  long,  but  the  slit  length  is  restricted  by  astigmatism. 

b)  In  the  case  of  Czerny -Turner  type  in  non-concentric  setting,  the  form  of 
isochromat  is  elliptic.  Therefore,  only  the  elliptic-an --shaped  slits  can  remove  the 
wavelength  error.  However,  when  the  coordinates  of  the  centers  of  the  spherical 
mirror  change,  the  form  of  the  ellipse  also  changes  rendering  the  adjustment  difficult. 

c)  In  Pfund  type,  long  slits  are  not  used.  Therefore,  if  there  is  no  loss  of  light 
energy  in  bringing  the  ray  onto  the  axis,  this  type  may  serve  well  in  the  near- 
infrared  region,  but  not  in  far-infrared  region. 

d)  Whichever  the  type  it  is,  the  effect  of  astigmatism  is  enhanced  as  the  slit 
becomes  long  which  rules  out  the  use  of  long  slits.  However,  this  effect  is  minimized 
by  setting  the  mirrors  and  the  grating  slightly  off-axis. 

e)  The  optimum  grating  size  is  proportional  to  and  approximately  to  /?t  in 
customary  setup. 
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Abstract 

A  vacuum  spectrophotometer  with  a  40  cm  concave  grating  has  been  con¬ 
structed  for  studies  of  solid  state  physics.  A  new  method  of  mounting  the 
grating  developed  by  one  of  the  authors  is  used  in  the  monochromator  in 
which  the  concave  grating  is  rotated  about  an  axis  somewhat  displaced  from 
the  center  of  the  grating.  With  this  arrangement,  the  resolution  is  found 
adequate  for  most  purposes  in  the  study  of  solid  state  physics.  This  spectro¬ 
photometer  can  be  used  for  transmittance  measurements  in  a  wide  wave¬ 
length  range  from  1050  A  to  6000  A  and  in  a  temperature  range  down  to  the 
liquid  nitrogen  temperature.  An  X-ray  irradiator  is  provided  for  coloration  of 
crystals.  Studies  of  fluorescence  and  phosphorescence  are  also  possible. 


1.  Introduction 

It  has  recently  come  into  notice  that  studies  of  optical  properties  of  solids  in  the 
vacuum  ultraviolet  region  made  important  contributions  to  the  advancement  of  solid 
state  physics.  The  spectrophotometers  used  in  those  previous  studies  were,  however, 
not  very  convenient  and  experimental  difficulties  sometimes  lowered  the  accuracy  of 
measurements  and  limited  the  field  of  investigation.  Development  of  spectrophoto- 
metric  techniques  is  at  present  one  of  the  most  important  requirements. 

A  vacuum  spectrophotometer  was  constructed  in  our  laboratory  and  several  devices 
used  were  proved  to  be  effective  in  photometric  work  on  color  centers  in  alkali  halides. 

2.  Monochromator 

High  resolving  power  of  monochromator  is  not  so  important  for  studies  of  solid 
state  physics  as  for  studies  of  gasious  molecules  because  of  the  inevitable  broadening 
of  spectral  lines  in  solid  materials.  In  most  cases  maximum  resolution  of  0.01  ev  is 
adequate,  which  corresponds  to  about  2  A  at  1500  A  and  3.5  A  at  2000  A  in  wavelength 
resolution.  For  a  vacuum  spectrophotometer,  smallness  of  volume  of  chamber  is  more 
desirable  than  high  dispersion  both  for  convenience  and  cost.  A  40  cm  concave  grat- 
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ing  with  600  grooves  per  mm  was  used  in  the  present  spectrophotometer,  and  the 
small  spectrophotometer  that  was  built  was  found  capable  in  separating  two  lines  3  A 
apart.  This  limit  of  resolving  power  is  not  due  to  the  smallness  in  size  of  the  monochro¬ 
mator,  but  to  insufficient  intensity  of  the  light  source  used. 

Several  types  of  grating  mounting  for  vacuum  monochromators  have  already  been 
developed.  As  to  the  mounting  condition  for  the  present  work,  the  following  are 
desirable  to  be  fulfilled. 

1.  Fixed  light  beam.  The  direction  of  the  light  beam  emerging  from  the  exit  slit  is 
to  be  fixed,  otherwise  the  light  beam  would  illuminate  different  points  on  the  speci¬ 
men  when  the  wavelengths  are  changed.  The  mounting  of  a  concave  grating  on  the 
end  of  an  arm  pivoted  at  the  center  of  the  Rowland  circle  used  by  Fujioka  and  Ito,*’ 
and  by  Tousey  et  al,'^  is,  therefore,  not  suitable  for  the  present  purpose  in  spite  of 
its  excellent  focusing.  Rotation  of  a  concave  grating  around  an  axis  at  or  near  the 
center  of  the  grating  is  desirable. 

2.  Good  efficiency.  Efficiency  in  energy  is  important,  for  the  reflective  power  of  the 
grating  surface  is  very  low  in  the  far  ultraviolet  region.  Except  the  case  of  extreme 
grazing  incidence,  normal  incidence  is  desirable,  for  the  astigmatism  is  small  and 
almost  full  energy  emerges  through  the  exit  slit. 

3.  Little  polarization  at  reflection  hy  the  grating.  For  this,  normal  incidence  is  again 
preferred. 

4.  Sharp  focusing.  Sharp  focusing  becomes  important  in  the  short  wavelength  range, 
for  the  required  resolution  is  necessarily  high.  Resolution  of  0.01  ev  corresponds  to 
0.8  A  at  1000  A,  but  to  20  A  at  5000  A.  In  the  long  wavelength  range  it  matters  little 
even  if  the  focusing  is  not  sharp. 

In  the  present  spectrophotometer,  the  concave  grating  is  rotated  about  a  vertical 
axis  displaced  a  little  from  the  center  of  the  grating.  Blurring  of  image  due  to  rotation 
is  reduced  by  simultaneous  translation  of  the  grating.  This  mounting  of  grating  was 
invented  by  Johnson®’  and  later  'improved  by  one  of  the  present  authors.^’  The 
monochromator  is  shown  schematically  in  Fig.  1.  The  entrance  and  exit  slits  Si,  Sj 
are  fixed  in  the  directions  making  an  angle  of  about  30  degrees  at  the  concave  grat¬ 
ing  G.  The  axis  of  rotation  of  the  grating  is  at  O,  displaced  on  the  plane  of  rotation 
perpendicularly  to  bisecter  of  the  angle  SiGSj.  The  optimum  displacement,  with  which 
the  monochromator  can  maintain  good  focus  as  the  grating  is  rotated,  is  given  by  the 

1.  Fujioka  and  Ito,  Science  of  Light,  1  (1951)  1. 

2.  Tousey,  Johnson,  Richardson  and  Toran,  J.  Opt.  Soc.  Am.,  41  (1951)  696. 

3.  P.  D.  Johnson,  Rev.  Sci.  Inst.,  28  (1957)  833. 

4.  Ryumyo  Onaka,  Science  of  Light.  7  (1958)  2.3, 
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following  equation 


/=/?sin^/{l— tan?i(tanao— tan^o)/2)  , 

where  / :  optimum  displacement  of  the  center  of  rotation, 

0 :  angle  between  the  grating  normal  and  the  bisecter  of  segment  Si  S,, 

^  ;  half  of  the  angle  SiGSj, 

Uo.  angle  of  incidence  monochromator  is  set  for  the  wavelength 

angle  of  diffraction 

taken  as  standard. 

These  quantities  are  related  one  another  as 


ao=0+ijt , 

$0  =  0  —  ^  . 

The  sine-bar  scanning  mechanism  A  (Fig.  1)  turns  the  shaft  which  is  directly  con¬ 
nected  through  a  vacuum-tight  sleeve  to  the  turn  table  for  the  grating.  The  wave¬ 
length  of  emerging  light  from  the  monochromator  is,  therefore,  linearly  proportional 


G :  Grating ;  O :  Center  of  Rotation ;  A :  Sine-bar  Scanning 
Mechanism;  C:  Wavelength  Counter;  M :  Synchronous  Motor 
S, :  Entrance  Slit ;  S^ :  Exit  Slit ;  L :  Light  Source ;  D :  Detector 
Head;  P:  Phototube  Housing;  X:  X-ray  Tube;  T„  T^,  T|: 
Cold  Traps  and  Diffusion  Pump. 
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to  the  advance  of  the  scanning  screw  and  is  indicated  on  a  digital  counter  in  Angstrom 
units. 

The  entrance  and  exit  slits  are  both  of  symmetrically  opening  type,  the  width 
of  which  are  made  adjustable  from  outside  of  the  vacuum  chamber.  Between  the 
slits  and  the  monochromator  body,  large  stop  cocks  are  provided  to  facilitate  disjoining 
of  the  light  source  or  the  detector  head  from  the  vacuum  chamber  whenever  needed. 

Three  vacuum  systems  are  incorporated,  one  each  for  the  monochromator,  the 
detector  head  and  the  additional  equipment  consisting  of  an  X-ray  tube,  an  evaporation 
chamber  etc.  The  vacuum  system  for  the  monochromator  comprises  a  liquid  nitrogen 
trap,  a  3  inch  oil  diffusion  pump  and  a  rotary  oil  pump. 

For  measurements  at  wavelengths  longer  than  1050  A,  a  lithium  fluoride  window 
is  provided  to  isolate  the  detector  head  without  impairing  its  high  vacuum.  This 
allows  the  monochromator  chamber  to  be  filled  with  pure  hydrogen  at  appropriate 
pressure  for  the  operation  of  the  light  source.  The  hydrogen  gas  seems  to  protect 
the  grating  against  sputtering  of  the  light  source. 

3.  Light  Source 

As  the  light  source  a  hydrogen  arc  discharge  lamp  was  used,  which  emits  sufficient 
energy  from  below  1000  A  to  visible  region.  The  quartz  capillary  tube,  6  mm  in  inner 
diameter  and  120  mm  long,  and  the  pure  aluminium  electrodes  of  the  lamp  are  all 
water-cooled.  0-rings  made  of  silicon  rubber  were  used  to  seal  the  joints  between  the 
quartz  capillary  and  the  electrodes.  The  light  source  is  operated  at  700  volt  and  about 
0.3  amp  dc  through  a  half  wave  rectifier  from  ac  high  voltage  supply.  This  gives  50 
cycle  ripples  to  the  light  output  which  are  beneficial  to  ac  detection  of  the  phototube 
response  enabling  the  discrimination  of  transmitted  light  from  the  long  life  phos¬ 
phorescence  of  crystal  samples. 

4.  Detector  Head 

The  detector  head  was  devised  for  the  following  measurements  and  treatments. 

1.  Transmittance  measurement. 

2.  Measurement  of  fluorescence  and  phosphorescence. 

3.  X-ray  coloration. 

4.  Photoconductivity  measurements. 

5.  Sample  making  by  vacuum  evaporation. 

(4  and  5  are  in  prospect.  The  present  X-ray  tube  is  intended  to  be  replaced  by  an 
evaporation  equipment  when  needed.) 

The  detector  head  is  shown  in  Fig.  2.  For  transmittance  measurements,  the 
specimen  crystal  held  in  a  holder  is  placed  in  contact  with  the  bottom  of  a  cryostat 
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and  can  be  cooled  down  to  the  liquid  nitrogen  temperature.  The  position  of  the  cry¬ 
stal  is  about  7  mm  away  from  the  axis  of  the  cryostat  cylinder,  which  allows  the 
crystal  to  be  in  or  out  of  the  light  beam  by  rotating  the  cryostat  The  cryostat  is 
sealed  vacuum  tight  by  means  of  an  0-ring.  Atmospheric  pressure  on  the  cryostat  is 
borne  by  a  thrust  bearing. 

Two  positions  are  possible  for  the  crystal  to  be  in  the  light  beam,  near  to  and  far 
from  the  light  detector.  The  crystal  faces  the  light  detector  (fluorescent  screen)  with 
different  solid  angle  at  these  two  positions,  and  scattered  light  at  the  crystal  surface 
and  fluorescence  of  the  crystal,  if  any,  give,  therefore,  two  different  contributions  to 
the  response  of  the  detector.  The  radio  between  these  two  components  in  the  response 
depends  only  upon  the  geometrical  relation  between  the  detector  and  the  crystal.  Thus 
it  is  possible  to  find  this  excess  of  response  from  two  readings  of  the  detector.  The 
ratio  of  the  excess  component  by  one  reading  to  that  by  the  other  was  approximately 
one  to  four  both  in  theoretical  estimation  and  experimental  determination  by  the  use 
of  a  fluorescent  crystal.  The  true  transmittance  is  therefore  found  if  a  third  of  the 
difference  between  the  two  readings  obtained  at  above  two  positions  is  substracted 
from  the  transmittance  obtained  at  the  farther  position  from  the  detector. 

A  combination  of  a  fluorescent  screen  and  RCA  IP 28  photomultiplier  was  used 
as  the  detector.  As  the  fluorescent  material,  sodium  salicylate  was  used,  and  the 


X  -  rajr  $ 


Fig.  2.  Detector  Head.  S:  Exit  Slit,  Wi :  Li  F  Window,  CR:  Cry- 
'  stal  Specimen,  F, :  KBr-Filter,  F,:  Glass- Filter,  FL:  Fluores¬ 
cent  Material,  W*:  Quartz  Window,  PH :  Photomultiplier  1P28 
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quartz  window  of  the  photomultiplier  housing  was  coated  with  it. 

The  output  of  the  photomultipler  is  fed  to  a  dc  or  ac  amplifier  and  then  to  an 
electronic  recorder.  When  the  sample  crystals  had  phosphorescence  as  the  result  of 
X-ray  irradiation  for  example,  the  ac  amplifier  tuned  to  50  cycle  was  used.  Since  ^e 
light  source  had  50  cycle  ripples,  50  cycle  component  of  the  transmitted  light  enabled 
the  true  transmittance  of  the  specimen  to  be  obtained,  for  the  afterglow  becomes 
undetectable  when  an  ac  amplifier  is  used. 

When  the  scattered  light  in  the  monochromator  is  negligible  as  in  the  case  of 
wavelengths  longer  than  1350  A,  addition  of  a  logarithmic  amplifier  facilitates  the 
calculation  of  transmittance. 

With  a  2  inch  oil  diffusion  pump  and  a  liquid  nitrogen  trap,  pressure  of  the  order 
of  10’^  mm  Hg  is  attained  in  the  detector  head. 

When  the  vacuum  spectrophotometer  is  used  for  near  ultraviolet  or  visible  region, 
use  of  a  set  of  filters  becomes  necessary  to  separate  higher  order  spectra.  The  follow¬ 
ing  filters  were  used. 

Up  to  1050  A  No  filter  No  window 

from  1050  A  to  2100  A  No  filter  LiF  window 

from  2100  A  to  3500  A  KBr  filter  LiF  window 

from  3500  A  to  6000  A  Glass  filter  LiF  window 

These  filters  were  mounted  on  a  cylindrical  holder  and 
were  made  exchangeable  by  rotating  the  holder  shaft 
from  outside  of  the  detector  head. 

4.  Fluorescence  and  Phosphorescence  Measurement 

Fluorescence  and  phosphorescence  of  crystal  samples 
can  be  measured  if  the  fluorescent  screen  on  the  photo¬ 
multiplier  housing  is  replaced  by  a  clean  quartz  window. 

The  spectra  of  the  fluorescene  can  be  taken  by  a  small 
quartz  spectrograph  with  a  mirror-lens  adaptor  instead 
of  the  photomultiplier  as  shown  in  Fig.  3. 

5.  X-Ray  Irradiator 


An  unsealed  X-ray  tube  with  a  tungsten  target  was 
provided  on  the  side  of  the  detector  head  for  studies  of 
color  centers  in  alkali  halides.  A  Hickman  oil  diffusion 
pump  was  used  for  evacuation.  The  X-ray  tube  is 
operated  at  30  KV  25  mA  in  most  cases. 


Fig.  3.  Lens-Mirror  System 
for  Studies  of  Luminescenc. 
CR:  Crystal  Specimen 
M :  Surface  Mirror 
L:  Quartz  Lens 
S:  Entrance  Slit  of  Spec¬ 
trograph. 
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The  spectrophotometer  can  be  used  over  a  wide  wavelength  range  from  1050  A  to 
6000  A  with  a  lithium  fluoride  window  and  a  hydrogen  discharge  tube.  Wavelength 
accuracy  was  within  ±  2  A  in  error  throughout  the  above  wavelength  range,  and  the 
best  resolution  obtained  was  about  3  A  at  the  vacuum  region  (below  20(X)  A). 


Fig.  5.  Transmittance  measurement  of  X-ray  irradiated  KBr  with  the  logarithmic  amplifier. 

A  typical  record  of  the  hydrogen  light  source  is  shown  in  Fig.  4,  and  that  of 
transmittance  measurement  with  the  logarithmic  amplifier  is  shown  in  Fig.  5. 

The  spectrophotometer  has  been  in  use  for  about  one  year  for  studies  of  color 
centers  in  alkali  halides  and  of  silver  activated  phosphors  of  NaCl.  Reports  on  these 
special  subjects  will  be  published  in  separate  papers.*^ 
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Summary 

In  the  previous  paper  it  was  stated  that  the  color  of  natural  object  is  one  of  the 
important  factors  in  treating  the  problem  of  visual  range  in  the  daytime,  and  the 
method  and  results  of  measurement  were  reported.  In  this  paper  the  visual  range 
of  colored  point  sources  is  discussed. 

From  the  comparison  of  measured  and  calculated  data  about  the  change  in 
apparent  luminance  and  chromaticity  of  natural  objects  with  distance,  the  extinction 
coefficient  of  light  in  the  atmosphere  is  estimated.  Extinction  coefficient,  ax,  is  ex¬ 
pressed  by  aoiioli)",  where  io  is  taken  to  be  0.55 Then  the  apparent  luminance  of 
natural  objects  is  calculated  by  replacing  the  mean  effective  extinction  coefficient  by 
<To  with  a  precision  of  two  percent  in  error.  Examples  of  the  estimation  show  that 
aA=2.5(io/x)®-*km~‘  on  a  rainy  day  when  the  visual  range  was  about  15  km,  and 
ffA=0.04(^il)*-^km"‘  on  a  clear  day. 

A  graphical  method  to  determine  the  visual  range  of  colored  light  sources  in  the 
daytime  is  developed.  The  results  gained  by  this  method  show  that  red  light  and 
green  light  of  adequate  saturation  have  good  visual  ranges  but  a  pale  blue  light  has 
poor  visibility  when  it  is  viewed  against  the  background  of  natural  objects. 


1.  Introduction 

In  discussing  the  visual  range  of  colored  light  sources  in  the  daytime,  important 
factors  to  be  considered  are  the  recognition  threshold  of  the  illuminance  produced  by 
the  sources  at  the  observer’s  eye,  the  luminance  and  color  of  background,  therefore 
the  colors  of  natural  objects,  and  the  extinction  coefficient  of  light  in  the  atmosphere. 
About  the  threshold  illuminance  of  colored  sources  the  author  already  reported*’,  and 
the  colors  of  natural  objects  were  dealt  with  in  the  previous  work**. 

1)  H.  Masaki :  J.  Ilium.  Eng.  Inst.  (Tokyo).  39,  (1955)  430. 

2)  II.  Masaki :  Science  of  Light  (Tokyo).  8,  (1959)  67. 
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The  extinction  coefficient  is  the  most  difficult  factor  to  obtain  experimentally. 
It  is  known  that  the  change  in  apparent  colors  of  natural  objects  with  distance 
relates  to  <Ji.  If  this  relation  can  be  used  for  the  estimation  of  at,  it  will  be  a  great 
advantage  in  treating  the  problem  of  visibility. 

A  graphical  method  to  estimate  a\  by  simply  measuring  the  apparent  colors  of 
objects  at  various  distances  was  introduced  and  applied  in  computing  the  visual  range. 

2.  Change  in  colors  of  natural  objects  with  distance 

The  apparent  colors  of  natural  objects  change  from  yellow-green  or  yellow  to  blue 
with  increasing  distance  and  finally  approach  the  color  of  the  horizon  sky.  The 
luminances  also  approach  that  of  the  sky.  This  is  because  the  light  from  objects 
attenuates  in  the  atmosphere  and  the  daylight  scattered  in  the  path  is  added  to  it. 
The  equation  expressing  the  apparent  luminance  Br  of  the  object  observed  at  the 
distance  r  was  first  derived  by  Koschmieder  and  arranged  to  be  more  comprehensive 
by  Duntley*’  and  Middleton^-  It  is 

Br  =  Bk{X-  e-"”)  +  Btfi-”',  .( 1 ) 

where  is  the  luminance  of  the  horizon  sky,  Bo  the  luminance  of  the  object  observed 
at  close  distance,  and  a  the  extinction  coefficient.  The  first  term  of  the  right  side 
expresses  the  contiibution  of  the  scattered  daylight  to  the  luminance  and  the  second 
term  the  contribution  of  the  attenuated  light  from  the  object.  Applying  the  equation 
(1)  to  a  monochromatic  light  of  wavelength  x,  we  obtain  the  following  equation 

Br,A  =  Bk,k{1  -  e->’^^)^Bo,ie-’'^^ .  (2) 

Middleton  calculated  the  change  of  colors  by  integrating  the  equation  (2)  over 
the  wavelength.  He  assumed  that  the  light  from  the  horizon  sky  had  the  same 
spectral  distribution  as  the  prevailing  illumination  and  expressed  Bo,  a  in  the  form 

Bo,  A  =  /a  (Cip  -f-  1)  Bh,  A  ( 3 ) 

where  /a  is  the  spectral  reflectance  of  the  object  and  Cw  the  contrast  between  tlie 
luminance  of  an  ideal  white  object  in  the  vicinity  of  the  observer  and  that  of  the 
horizon  sky.  As  for  the  spectral  intensity  distribution  of  the  sky,  that  of  the  CIE 
source  C  was  used,  which  represents  well  the  general  atmospheric  conditions.  But 
when  the  prevailing  illumination  and  the  light  from  the  horizon  sky  have  different 
spectral  qualities,  for  instance  when  on  a  clear  day  the  object  is  illuminated  by  the 

."?)  S.  Q.  Duntley :  J.  Opt.  Soc.  Amer.  38,  (1948)  179. 

1)  W.K.K.  Middleton :  J.  Opt.  Soc.  Amer.  40.  (1950)  373. 

5)  W.E.K.  Middleton:  Vision  Through  the  Atmosphere,  Toronto  Univ.  Press,  1952. 


Apparent  Colors  of  Natural  Objects  (//) 


41 


direct  sunlight  of  the  color  temperature  5  500“K  and  the  light  from  the  horizon  sky 
has  the  color  temperature  of  more  than  10,000’K,  the  following  more  precise  ap¬ 
proximation  is  obtained  by  integrating  (2)  directly : 

A'r=|  {Ba,a(1— ,  etc.  (4) 

In  Fig.  1  the  changes  in  the  chromaticity  of  foliage  (Xo=0.3292,  yo=0.3778)  against 
the  blue  sky  (xa =0.2750,  |/*=0.2954)  and  the  sky  of  the  same  spectral  intensity 


Fig.  1.  Spectral  distributions  and  chromaticity  coordinates  of 
the  light  from  foliage  at  various  distances  in  non-selective 
(fuli  lines)  and  moderately  selective  (dashed  lines)  atmos¬ 
pheres  (calculated  results). 

distribution  as  the  source  C,  calculated  by  (4)  for  non-selective  and  moderately  selec¬ 
tive  air,  are  shown  The  change  in  the  spectral  distribution  by  distance  is  also 
shown.  It  might  not  be  reasonable  to  consider  the  sky  to  be  blue  despite  the  atmos¬ 
phere  being  non-selective,  but  it  is  only  shown  as  an  example  of  calculation.  The 
calculated  results  show  that  the  change  is  greater  against  the  blue  sky  than  against 
the  sky  with  the  characteristic  of  the  source  C.  It  is  worth  noting  that  the  original 
color  of  the  object  is  lost  at  a  small  distance  even  in  a  clear  atmosphere. 

Use  of  an  electronic  computer  saved  the  otherwise  laborious  work  of  numerical 
integration  of  (4). 
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3.  Estimation  of  extinction  coefficient 

The  general  method  to  obtain  the  extinction  coefficient  a  of  light  in  the  atmos¬ 
phere  is  first  to  measure  the  transmittance  r,  through  the  layer  of  air  of  the 
thickness  r,  then  compute  by  the  equation 

Tr  =  .  (  5  ) 

By  measuring  the  transmittance  spectrophotometrically,  the  extinction  coefficient  ax 
of  light  for  the  wavelength  i  is  obtained*'^) .  But  this  method  is  attended  with 
several  difficulties.  First,  it  needs  an  experimental  arrangement  on  a  large  scale. 
Second,  the  air  is  often  unsteady;  variation  of  the  transmittance  during  scanning 
through  the  spectrum  causes  errors.  To  measure  the  transmittance  of  fog,  the  author 
designed  a  spectrotelephotometer  with  which  the  scanning  of  the  visible  spectrum  is 
performed  in  1/20  second,  the  spectral  transmittance  curve  being  shown  on  a  cathode 
ray  oscillograph.  But  this  instrument  cannot  be  used  over  a  long  distance  necessary 
to  measure  the  air  of  small  extinction  coefficient.  The  third  difficulty  is  that  most 
of  the  radiation  detectors  have  low  sensitivity  in  the  ends  of  visible  spectrum,  the 
violet  and  red  regions,  and  the  spectral  selectivity  of  ax  cannot  be  determined 
precisely. 

When  the  spectral  energy  distribution  Br,\  of  the  light  from  objects  at  various 
distances  is  measured,  ax  is  calculated  from  (2)  by 

ax={llr\ogioe)logio[iBx,x—Bo,x)KBK,i-Br,i)\.  (6) 

This  method  is  also  encountered  with  the  same  difficulties  as  before. 

If  the  extinction  coefficient  and  its  dependence  on  wavelength  are  estimated  from 
the  measurement  of  only  the  apparent  luminance  and  chromaticity  of  objects  at  various 
distances,  it  becomes  possible  to  obtain  them  for  a  variety  of  atmospheres  because  of 
the  experimental  simplicity.  For  this  purpose,  equations  (4)  are  converted  to  the 
summation  in  the  form 

Xr=E{^in+{Boi—Bki)e-''*^]ii,  etc.  (7) 

i 

When  Bxi  and  Boi  are  known  and  Xr,  Yr  and  Zr  are  measured  at  more  than  i/S  dis¬ 
tances,  <T<  is  determined  by  solving  the  simultaneous  equations  (7)  numbering  i  which 
is  unknown  at  present.  It  is  not,  however,  easy  to  solve  such  equations  even  by  an 
electronic  computor.  Besides,  the  measured  values  of  Xr,  Yr  and  Zr  may  contain 
experimental  errors,  hence  the  mean  values  of  are  to  be  determined  by  numbers 
of  measurements  as  many  as  several  times  of  i/S. 

Fortunately  it  is  known®’  that  ax  or  at  is  approximately  inversely  proportional  to 

6)  A.  Arnulf,  J.  Bricard,  E.  Cur^  et  C.  Veret :  Rev.  Opt.  38  (1959)  103. 

7)  H.  Masaki:  Railw.  Tech.  Res.  Inst.  Interm.  Rep.  7-120,  (1956). 
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a  power  of  the  wavelength  i  and  the  exponent  n  lies  in  the  range  of  about  0—4. 
at  is  then  expressed  as 

Oi  =  </o  •  (  8  ) 

Applying  (8)  to  (7)  we  obtain 


Bri= 

Xr=^BriXi,  Yr=2^BriVit  Zr=J^BrlZi 
I  i  i 


(9) 


where  the  units  of  Bm  and  Bri  are  so  determined  that  Yr  expresses  the  luminance 
of  the  object  in  cd/m*.  For  >io  the  wavelength  0.55  u  is  taken  because  it  is  nearly  the 
center  of  the  visible  spectrum  and  055  is  a  convenient  number  for  computation.  In 
(9),  the  unknown  number  i  in  (7)  reduces  to  two.  Thus  the  number  of  measure¬ 
ments  necessary  to  obtain  mean  values  is  much  reduced. 

When  the  mean  effective  extinction  coefficient  <;«  is  defined  here  by  the  following 
equation 

Br=Bk{l-e-<'»')+Boe-<'’^'  (10) 


Table  1.  Apparent  luminance  of  foliage  at  various 

distances  (calculated  results).  unit:  cd/m* 


(km)  ,^| 

0 

1 

2 

3 

4 

Case  1, 

=0.040  km*' 

0 

1150 

1.3 

1349 

1.350 

1352 

1.356 

1361 

11 

2545 

2551 

2560 

2576 

2589 

29 

3839 

3844 

3847 

3851 

3846 

oo 

5050 

Case  11.  tfu=0.045km*' 

0 

420 

0.3 

.559 

559 

560 

561 

563 

5 

2527 

2522 

2525 

2542 

2556  • 

21 

6814 

6796 

6778 

6772 

6745 

36 

8809 

8787 

8752 

8717 

8653 

OO 

10890 

Case  III, 

<»o=0.200  km*' 

0 

240 

0.3 

353 

353 

353 

354 

355 

0.75 

512 

511 

511 

513 

515 

3.5 

1222 

1218 

1215 

1216 

1213 

5.0 

1473 

1468 

1464 

1461 

1454 

OO 

2190 

Case  IV, 

tfj=2.5  km*' 

0 

130 

0.05 

198.2 

197.8 

197.8 

198.3 

198.7 

0.3 

436.3 

4.34.9 

433.6 

433.3 

432.0 

0.75 

621.6 

619.8 

617.3 

614.5 

610.0 

1 

663.0 

661.5 

659.0 

656.1 

651.5 

OO 

711 
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it  is  found  from  the  result  of  more  than  hundred  calculations,  some  of  which  are 
shown  later  (Table  1),  that  the  error  in  calculated  values  of  Bt  caused  by  replacing 
Om  by  <7o  in  (8),  the  extinction  coefficient  for  the  wavelength  0.55/1,  is  less  than  two 
percent.  This  happens  because  the  spectral  distribution  of  the  light  from  natural 
objects  is  not  conspicuously  selective.  By  putting  is  obtained  only  from 

the  change  in  luminance  with  the  distance,  namely 

ffo=(l/r  logio  e)  logio  {(BM-Bn)/(BM-Br))  (11) 

It  is  more  easily  estimated  graphically  as  shown  in  Examples. 

The  exponent  «  in  (8)  represents  the  degree  of  dependence  of  <7x  or  on  wave¬ 
length.  As  it  is  considered  to  be  about  0  —  4  in  the  ordinary  atmosphere,  chromaticity 
coordinates  of  the  light  from  objects  at  various  distances  are  calculated  by  (9)  with 
<7o  determined  by  (11)  and  «=0,  1,  2.  3,  and  4.  It  needed  somewhat  laborious  work 
to  calculate  manually.  But  once  a  program  for  the  electronic  computor  was  for¬ 
mulated,  many  computations  could  be  made  in  a  short  time.  Then  by  comparing  the 
measured  data  with  the  calculated  on  the  chromaticity  diagram.  «  is  easily  deter¬ 
mined.  Examples  are  shown  in  Table  1  and  Figs  2~5. 
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Examples 

Case  I.  On  a  clear  day  in  winter,  the  horizon  sky  looked  blue.  From  Fig.  2 
<To=0-04km"*  and  «=2.7  are  estimated. 

Case  II.  On  a  clear  day  in  spring,  the  horizon  sky  looked  white.  From  Fig.  3 
<7o=0.045  km'*  and  n=0  are  estimated. 

Case  III.  On  an  overcast  and  hazy  day  in  spring,  the  visual  range  was  about 
10  km.  From  Fig.  4  <To=0.20km'*  and  w=1.3  are  estimated. 

Case  IV.  On  a  rainy  day  in  spring,  the  visual  range  was  about  1.5  km.  From 
Fig.  5  ffo=2.5  km'*  and  m=0.5  are  estimated. 


Fig.  5.  Luminances  and  chromaticities  of  foliage  at  various  distances  (case  IV). 
Symbols  with  heavy  outlines  represent  experimental  results. 


4.  Visual  range  of  colored  point  sources  in  the  daytime 

The  illuminance  produced  by  a  light  source  at  the  observer’s  eye  decreases  in 
proportion  to  the  inverse  square  of  the  distance  and  by  the  extinction  through  the 
atmosphere.  The  maximum  distance  at  which  an  observer  can  recognize  the  color 
of  the  source  is  the  distance  at  which  the  illuminance  decreases  to  its  chromatic 
threshold  value. 
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The  illuminance  E  for  the  observer  at  the  distance  r  from  the  source  of  the 
intensity  /o  is  given  by 

E=Iae-’"lr^-=Irlr*  (12) 

where  E  is  the  apparent  intensity  of  the  source  observed  at  the  distance  r  if  the 
effect  of  the  scattered  daylight  is  neglected.  As  the  colored  sources  have  generally 
more  selective  spectral  distributions  than  the  light  from  natural  objects,  the  mean  ex¬ 
tinction  coefficient  a  cannot  be  replaced  by  oo  of  the  wavelength  0.55  /t.  It  is  calculated 
by 

a=  -(1/r  log,o«)  logio|  j  j  =  -(1/r  logio^)  logio(/r//o)  (13) 

where  Ei  is  the  spectral  intensity  distribution  of  the  source. 

The  chromatic  threshold  illuminance  Ete  is  expressed  as  pEt,  where  Et  is  the 
achromatic  threshold  value.  The  ratio  of  chromatic  to  achromatic  threshold  illuminance, 
p,  is  related  to  the  background  luminance  and  the  chromaticity  difference  between 
source  and  background*\  Expressing  the  chromaticity  difference  by  the  distance  on 
the  Breckenridge  and  Schaub’s  rectangular  uniform  chromaticity  scale  (RUCS)  diagram"’ 
the  author  determined  the  relation  as 

</(/»-!)*  =  a  (14) 

where  d,  the  distance  on  the  RUCS  diagram,  is  computed  by 

rf={(x/-x.')*+(y/-p*-')»)^  (15) 

from  the  RUCS  coordinates  of  the  source  (x/,  y,")  and  the  background  (xd",  Vb"). 
The  constant  a  is  0.101  for  the  group  of  sources  that  give  red  or  yellow  response  to 
the  observer  (red-yellow  recognition)  and  0.051  for  the  group  of  sources  that  give 
green  or  blue  response  (green-blue  recognition)  for  50%  detection  in  the  rangd  of 
background  luminance  1,0(X)— 10,0(X)cd/m*  corresponding  to  the  luminance  in  the 
daytime.  Et  is  the  achromatic  threshold  value  of  red  light  for  red-yellow  recognition 
or  of  green  light  for  green-blue  recognition  against  an  achromatic  background.  In 
the  same  range  of  background  luminance,  the  achromatic  threshold  value  in  Im/km* 
for  50%  detection  is  given  by 

logic  £» = 0.9  logic  Bb—b  (16) 

where  Bb  is  the  background  luminance  expressed  in  cd/m*,  and  b  is  1.30  for  red  light 
and  1.07  for  green  light.  The  relations  Et  vs.  Bb  and  p  vs.  d  are  shown  in  Fig.  6. 
In  open  air  observation,  natural  objects  make  the  background ;  mean  values  of 


8)  F.  C.  Breckenridge  and  W,  R.  Schaub :  J.  Opt.  Soc.  Amer.  29  (1939)  370. 
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Fig.  6.  Achromatic  threshold  illuminances  of  red  and  green  point  sources.  Et, 
as  the  funtion  of  background  luminance.  Bb,  and  ratios  of  chromatic  to 
achromatic  thresholds  for  red-yellow  and  green-blue  recognitions  against 
background  of  luminance  1000— 10000  cd/m*,  p,  as  the  function  of  chroma- 
ticity  difference  between  source  and  background,  d. 


Table  2.  Mean  values  of  apparent  luminance  of  natural  objects  (measured  data). 


Object 

Distance 

(km) 

Season 

Weather 

Luminance 

(cd/m2) 

foliage 

0-1  1 

spring 

clear 

530 

overcast 

240 

summer 

clear 

710 

overcast 

380 

i 

autumn 

clear 

550 

overcast 

300 

winter 

clear 

730 

overcast 

550 

6-«  -I 

spring 

clear 

1800 

'  ! 

overcast 

450 

16  i 

spring 

clear 

2050 

grass 

0.1  1 

summer 

overcast 

880 

winter 

clear 

1580 

field  crops 

0.1  i 

summer 

clear 

1340 

overcast 

370 

earth 

0.1 

winter 

clear 

1100 

overcast 

510 

sea  surface 

1  j 

winter 

clear 

1840 

1 

overcast 

880 

snow  surface 

0.1  ! 

winter 

clear 

12400 

building 

2 

spring 

overcast 

470 

horizon  sky 

winter 

clear 

7100 

spring 

overcast 

2000 
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their  luminance  are  shown  in  Table  2.  The  average  values  of  their  chromaticity 
coordinates  reported  in  the  previous  paper  are  transformed  to  the  RUCS  coordinates 
and  shown  in  Table  2. 

Equating  E  given  by  (12)  to  Ete  or  pEt,  we  obtain  the  visual  range  r,  from 
the  equation 

E,lh=e-<''-lprr*  (17) 

For  any  combination  of  the  source  and  background,  e'^'lpr*  is  computed  by  (13), 
(14)  and  (15)  and  the  e'^'^jpr^  vs.  r  curve  is  plotted  on  a  diagram.  Given  the  values 
of  Bb  and  To,  Etlh  can  be  computed  from  (16)  and  the  visual  range  can  be  obtained 
from  the  curve. 

As  examples,  four  sources,  red  (S 1),  yellow  (S  2),  green  (S  3)  and  pale  blue 
(S4)  are  used.  Their  spectral  intensity  distributions  are  shown  in  Fig.  8.  As 
the  state  of  the  atmosphere  three  cases  are  given :  case  V  (<to=0.02  km'*,  «=3),  case 
IT  (<to=0.20  km'*,  «=2)  and  case  IIT  («To=2.0km"*,  n=l).  The  backgrounds  are 


Fig.  7.  Mean  values  of  chromaticities  of  natural 
objects  (RUCS  coordinates). 


WAVELENGTH  (mM) 


Fig.  8.  Spectral  intensity  distribution  curves  of  red  (S  1) 
yellow  (S2),  green  (S3),  and  pale  blue  (S4)  lights. 


Fig.  9.  RUCS  coordinates  of  lights  attenuated  in  the 
atmosphere  (calculated  results). 
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Table  3.  Ratio  of  the  apparent  intensity  of  light  sources 
observed  at  distance  r  to  the  initial  intensity, 

/r//o<  (calculated  results). 


n 

0 

1 

2 

3 

"1 

0.01 

S  1  (red) 

0.990 

0.992 

0.993 

0.993 

0.10 

0.905 

0.919 

0.928 

0.939 

0.50 

0.607 

0.650 

0.689 

0.724 

1.00 

0.369 

0.422 

0.474 

0.524 

2.50 

0.082 

0.116 

0.154 

0.200 

5.00 

0.007 

0.013 

0.024 

0.042 

0.01 

S2  (yellow) 

0.990 

0.991 

0.992 

0.992 

0.10 

0.905 

0.909 

0.912 

0.915 

0.50 

0.609 

0.622 

0.635 

0.645 

1.00 

0.368 

0.385 

0.401 

0.416 

2.50 

0.083 

0.093 

0.105 

0.118 

5.00 

0.007 

0.009 

0.011 

0.016 

0.01 

S3  (green) 

0.991 

0.991 

0.990 

0.990 

0.10 

0.905 

0.904 

0.902 

0.900 

0.50 

0.608 

0.602 

0.596 

0.590 

1.00 

0.368 

0.363 

0.357 

0.350 

2.50 

0.082 

0.080 

0.079 

0.078 

5.00 

0.006 

0.006 

0.006 

0.006 

0.01 

S4  (pale  blue) 

-  0.992  0.992 

0.992 

0.991 

0.10 

0.906 

0.906 

0.905 

0.904 

0.50 

0.607 

0.607 

0.606 

0.604 

1.00 

0.369 

0.369 

0.369 

0.368 

2.50 

0.082 

0.083 

0.087 

0.090 

5.00 

0.007 

0.007 

0.009 

0.010 

foliage  in  the  vicinity  FI  (x= 0.329,  0.369)  and  a  distant  object  F2  (x= 0.282, 

If =0.301).  The  variation  of  the  apparent  intensity  and  chromaticity  of  the  sources 
in  the  atmosphere  are  shown  in  Table  3  and  Fig.  9.  The  e'^'^lpr'  vs.  r  curves 
are  shown  in  Fig.  10. 

With  red  light  fS  1),  its  chromaticity  is  very  different  from  that  of  the  back¬ 
ground  and  therefore  p  is  small.  Its  apparent  color  changes  very  little  by  the 
extinction  and  the  visual  range  is  good  in  various  states  of  the  atmosphere  and 
background.  Yellow  light  (S2),  though  not  much  affected  by  the  atmosphere, 
has  a  smaller  visual  range  against  the  background  of  foliage  or  earth  in  com¬ 
parison  with  red  light  of  the  same  intensity.  As  for  green  light,  it  is  fairly 
distinct  in  color  from  the  background.  Though  its  color  changes  by  the  selective 
extinction,  it  is  easily  recognized  because  the  locus  of  chromaticity  deviates  fairly 
from  natural  objects.  Green  light  of  saturation  lower  than  S  3  have  poor  visibility, 
for  they  become  assimilated  in  color  with  the  background.  Pale  blue  light  (S4) 
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DISTANCE  r  (km) 

Fig.  10.  e-''’^lpr^  as  the  function  of  distance  r.  Full  lines  represent  observa¬ 
tions  against  background  F  1  and  dashed  lines  observation  against  F  2. 

has  the  least  visual  range,  for  it  resembles  the  color  of  natural  objects ;  it  changes 
to  white  and  yellow  by  selective  extinction. 

Given  5o=2,500  cd/m*  and  /o=l,0(X)  cd,  £f//o=0.077  for  red-yellow  recogni¬ 
tion  and  0.133  for  green-blue  recognition.  Against  these  two  backgrounds,  F 1 
and  F2,  the  visual  ranges  obtained  from  Fig.  10  are  1.2-3.4km,  l.l-2.8km  and 
0.9-2.1km  for  SI,  S2  and  ^3  respectively.  For  S4  the  ranges  are  0.6-1.7km 
against  F 1  and  only  0.13-0.42  km  against  F  2. 

5.  Discussion 

In  the  estimation  of  extinction  coefficients  from  the  change  of  colors,  it  is  assumed 
that  the  spectral  reflectances  of  measured  objects  are  equal  and  the  atmosphere  is 
uniform.  If  equally  finished  targets  are  located  at  various  distances,  the  first  assump¬ 
tion  is  justified  imless  they  fade  differently.  But  such  arrangement  on  a  large  scale 
would  deprive  the  instrumental  simplicity  of  its  advantage.  In  the  examples  in  section 
3,  foliage  at  various  distances  was  actually  measured.  To  investigate  the  error  in¬ 
volved  in  this  method,  the  change  of  colors  is  computed  on  two  sets  of  initial  values. 
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Fig.  11.  Chromaticities  of  the  light  from  foliage  at 
various  distances  showing  the  effect  of  different 
initial  values  (calculated  results). 


Boi,  with  the  same  Bm  and  The  difference  between  the  initial  values  of  these 
two  sets  is  thought  to  be  maximum  for  foliage  of  the  same  kind  measured  at  the 
same  time.  The  results  are  shown  in  Fig.  11.  The  difference  is  fairly  large  at  a 
small  distance  but  decreases  rapidly  with  the  increase  in  distance.  It  is  found  that 
the  colors  of  distant  objects  are  affected  by  the  color  of  horizon  sky  and  the  charac¬ 
teristic  of  extinction  coefficients,  especially  their  dependence  on  wavelength  rather 
than  by  initial  colors.  As  the  computation  of  the  visual  range  requires  only  rough 
values  of  ai,  it  suffices  to  select  objects  of  the  same  kind  at  various  distances. 

In  the  computation  by  (9),  the  spectral  energy  distributions  of  light  from  the 
horizon  sky,  Bm,  and  from  object  near  the  observer,  B^u  are  necessary.  However, 
because  the  spectral  intensity  distributions  of  natural  objects  do  not  differ  with  one 
another  greatly  as  shown  in  the  previous  paper,  spectral  distributions  so  corrected  by 
the  already  known  distributions  for  the  object  of  the  same  kind  as  to  be  consistent 
with  the  measured  values  of  the  luminance  and  chromaticity  can  be  applied.  It 
causes  a  very  small  error.  Therefore  extinction  coefficients  are  estimated  by  the 
measurement  of  apparent  luminance  and  chromaticity  alone. 
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To  justify  the  second  assumption,  the  atmosphere  without  visible  layers  of  local 
cloud  or  fog  must  be  chosen. 

6.  Conclusion  and  acknowledgement 
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dence  on  wavelength  were  estimated.  By  applying  data  on  chromatic  threshold 
illuminance,  colors  of  natural  objects  dealt  with  in  the  previous  report  and  the  esti¬ 
mated  extinction  coefficient,  visual  range  of  colored  point  sources  in  the  daytime  was 
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Abstract 

The  infrared  spectra  of  HCl  associated  with  proton  accepters  in  the  region 
of  fundamental  stretching  vibration  reveal  two  broad  absorption  bands.  These 
bands  are  shown  to  arise  from  an  energy  level  splitting  and  can  be  interpreted 
satisfactorily  by  assuming  an  asymmetrical  double  potential  in  the  energy 
curve  for  the  proton.  The  result  of  observation  made  on  the  spectra  leads 
to  the  conclusion  that  the  second  minimum  lies  in  the  vicinity  of  unperturbed 
(p=l)  vibrational  level  of  the  first  minimum.  Further  the  new  absorption 
bands  observed  in  ultraviolet  region  relate  to  X-H-Cl  hydrogen  bond  and 
seem  to  be  attributed  to  the  charge  transfer  force. 


1.  Introduction 

As  regards  the  existence  of  a  second  minimum  in  potential  energy  curve  for  a 
proton  in  a  weak  hydrogen  bonded  system,  many  theoretical and  experimental*’ 
studies  have  been  carried  out  with  much  interest.  Sato^’  observed  two  new  bands  of 
the  first  OH  stretching  overtone  in  hydrogen-bonded  solvents  and  gave  a  detailed 
interpretation  by  assuming  a  double  minimum  potential  for  the  motion  of  proton. 
Recently  Barrow  and  his  coworkers*’  “’  observed  a  similar  splitting  of  the  first  OH 
stretching  overtone  of  alcohols  and  concluded  theoretically  that  a  second  potential 
minimum  lies  in  the  region  of  v=2  OH  vibrational  level.  However,  it  is  rather  a 

1)  C.  A.  Coulson  and  V.  Danielson,  Arkiv.  Fysik.  8  (1954)  239. 

2)  E.  R.  Lippincott  and  R.  Schroeder,  J.  Chem.  Phys.  23  (1955)  1099. 

3)  E.  R.  Lippincott,  J.  Chem.  Phys.  26  (1957)  1678. 

4)  C.  Reid.  J.  Chem.  Phys.  30  (1959)  182. 

5)  P.  C.  Mckinery  and  G.  M.  Barrow,  J.  Chem.  Phys.  31  (1959)  294. 

6)  Y.  Sato  and  S.  Nagakura,  Science  of  Light.  4  (1955)  120. 

7)  Y.  Sato,  Science  of  Light.  6  (1957)  108. 

8)  R.  Blinc  and  D.  Hadzi.  Mol.  Phys.  1  (1958)  391. 

9)  C.  L.  Bell  and  G.  M.  Barrow,  J.  Chem.  Phys.  31  (1959)  300. 

10)  C.  L.  Bell  and  G.  M.  Barrow,  J.  Chem.  Phys.  31  (1959)  1158. 
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question  whether  t>=0  or  v=\  is  the  region  in  which  the  second  potential  minimum 
occurs.  Although  thorough  examination  of  the  OH  fundamental  absorption  band  would 
clear  up  the  matter,  correct  measurment  in  this  region  is  difficult  because  of  the 
strong  intensity  of  CH  stretching  absorption  bands  in  solvents,  while  the  separation 
of  this  absorption  has  been  observed  in  OD  fundamental  band  of  phenols  which  is 
free  from  such  interfering  absorption. 

In  order  to  avoid  the  overlapping  of  CH  stretching  absorption  in  the  fundamental 
region  and  further  to  enlighten  us  on  the  question  of  the  height  of  the  second  poten¬ 
tial  minimum  in  X-H-Cl  system,  which  is  different  from  0-H-X  system  dealt  with  in 
the  former  studies,  an  experiment  was  made  on  the  HCl  fundamental  absorption  band 
in  proton  accepters.  The  remarkable  shift  of  HCl  vibrational  absorption  band  in  polar 
solvents  has  already  been  observed  and  interpreted  as  being  attributed  to  hydrogen 
bonding  by  Martin,  Gordy  and  Szobel,“’‘^®^  but  its  splitting  in  infrared  region  of  the 
spectrum  has  not  been  reported.  Taylor  and  Vidale*®’  found  two  absorption  bands  in 
the  Raman  spectrum  of  1 : 1  complex  of  HCl  and  dimethylether,  but  did  not  give  ap¬ 
propriate  interpretation.  While  measuring  HCl  fundamental  absorption  band  of  a  few 
proton  accepters  on  their  infrared  spectra,  the  auther  found  in  every  case  a  weak  second 
absorption  band  on  the  long  wavelength  side  of  the  strikingly  shifted  main  band. 
For  these  two  bands,  a  very  satisfactory  explanation  was  possible  by  assuming  an 
asymmetrical  double  potential  in  X-H-Q  hydrogen  bond  formed  between  molecules  of 
HCl  and  the  solvent.  The  second  minimum  was  found  to  lie  1800  cm"*~1900  cm“® 
higher  than  the  vibrationless  (t;=0)  level  and  lower  than  v=\  level  of  HCl  and 
capable  of  interacting  with  the  latter  level.  Between  the  height  of  this  minimum  and 
the  relative  basicity  of  proton  accepters,  a  characteristic  relation  was  found. 

In  describing  the  hydrogen  bonding,  the  ionic  and  covalent  nature  of  bonding 
electons  and  the  role  played  by  non-bonding  electrons  as  well  as  the  motion  of  protons 
should  obviously  be  taken  into  consideration.  On  the  absorption  spectum  of  the  complex 

t 

of  HCl  and  solvent  m  ultraviolet  region,  a  new  absorption  band,  which  is  not  observed 
either  in  HCl  or  in  the  solvent,  was  found.  The  charge  transfer  force  found  out  by 
Mulliken,“’"‘”  or  in  other  words,  the  transition  of  non- bonding  electrons  of  the  proton 
accepter  to  an  excited  state  level  of  HCl,  is  considered  responsible  for  this  new  ab¬ 
sorption  band. 

11)  W.  Gordy  and  P.C.  Martin.  J.  Chem.  Phys.  7  (1939)  99. 

12)  W.  Gordy  and  P.C.  Martin,  J.  Chem.  Phys.  9  (1941)  215. 

13)  L.  Szobel,  Comp.  rend.  218  (1944)  347,  843. 

14)  G.  L.  Vidale  and  R.  C.  Taylor,  J.  Am.  Chem.  Soc.  78  (1956)  294. 

15)  R.  S.  Mulliken,  J.  Chem.  Phys.  7  (1939)  20. 

16)  R.  S.  Mulliken,  J.  Am.  Chem.  Soc.  72  (1950)  600. 

17)  R.  S.  Mulliken,  J.  Am.  Chem.  Soc.  74  (1952)  811. 
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2.  Experimental  and  results 

For  obtaining  absorption  spectra  in  infrared  region,  a  Koken  I>-301  spectrometer 
with  two  NaCl  prism  and  25~50a*  rocksalt  cell,  and  in  ultraviolet  region,  a  Hitachi 
EPU-2  spectrometer  with  a  qurtz  prism  and  qurtz  cell  with  1  cm  path  length  were 
used.  In  infrared  region,  frequency  calibration  was  made  for  each  spectrum  by  super¬ 
posing  polystyrene  and  carbon  monoxide.  HCl  was  used  as  the  proton  donor  and 
several  compounds  such  as  ethers  and  ketones  with  lone  pair  electrons  as  the  proton 
accepters.  Molar  concentration  of  HCl  in  solvents  was  made  reasonably  low  for  1:1 
complex  to  be  formed.  Care  was  taken  of  impurities,  particularly  water  for  it  is  likely 
to  form  hydrogen  bond  with  the  solvent. 


Fig.  1.  1)  1.9  mol  HCl  in  diethylether 
2)  3.8  mol  HCl  in  dioxane 
Cell  thickness ;  25  ^ 


Table 


1  va  1 

( 

AV 

Diethylether 

2340  cm-‘  1 

1895  cm-* 

'  455  cm-' 

Isopropylether 

2370  1 

1966  ' 

404 

Dioxane 

2322 

1945 

377 

Acetone 

2319 

1969 

351 

Tetrahydrofuran 

2231 

1891 

331 
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On  infrared  absorption  spectra,  two  broad  bands,  one  strong  in  the  region  of 
2400  cm‘*  and  the  other  weak  in  the  region  of  1900  cm~S  were  always  observable. 
The  absorption  curves  of  HCl-dioxane  and  HCl-diethylether  solution  are  illustrated  in 
Fig.  1  as  examples,  and  the  wavelength,  separation  and  relative  intensity  of  the  two 
bands  appeared  in  other  cases  are  given  in  Table  with  curves  in  Fig.  4. 

Further,  on  dioxane  and  diethylether,  to  which  HCl  was  added,  a  sharp  band  was 
observed  at  837  cm”*  and  826  cm”*  respectively.  The  band  in  the  case  of  dioxane-HCl 
solution  is  shown  in  Fig.  2. 


Wavelength  (C7n”0 

Fig.  2.  1)  pure  dioxane 

2)  3.8  mol  HCl  in  dioxane 

3)  6.8  mol  HCl  in  dioxane 

Cell  thickness:  25 /i 

In  ultraviolet  absorption  spectra  of  proton  accepters  with  HCl,  a  broad  absorption 
band  was  found  in  the  neighborhood  of  360  m/^,  except  for  diethylether  and  isopropy- 
lether  of  which  the  maxima  were  obscure.  The  maximum  of  this  absorption  band 
was  at  365  m/t  for  methylethylketone,  at  a  little  above  360  m/i  for  acetone,  at  360  m/i 
for  ethylacetate,  at  a  little  below  360  m/i  for  dioxane.  On  the  other  hand  in  inactive 
solvent  carbon  tetrachloride,  the  absorption  band  was  at  the  wavelength  shorter  than 
300  mft  overlapping  the  characteristic  absorption  band  of  carbon  tetrachloride  as  clearly 
shown  in  Fig.  3. 
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Fig.  3.  a)  pure  CCI4  b)  a  small  amount  of  HCl  added  c)  pure 
acetone  d)  1.6  mol  HCl  in  acetone 
Cell  thickness ;  10  mm 


*  3.  Discussion 

The  cause  of  appearance  of  the  two  absorption  bands  in  2400  cm*  *  and  1900  cm*  ‘ 
regions  may  be  the  broadening  of  P  and  R  rotational  structures,  but  considering  the 
marked  shift  of  these  bands  from  2890  cm**  absorption  band  of  free  gasious  HCl,  we 
are  certain  that  the  hydrogen  bond  is  formed,  and  the  free  rotation  of  HCl  molecule 
becomes  therefore  not  possible  as  supported  by  various  other  evidence.  Should  how¬ 
ever  the  restricted  rotation  be  possible,  we  would  expect  a  linear  relation  between  the 
separation  of  the  bands  and  the  dielectric  constant  of  the  solvent.  The  dielectric  con¬ 
stant  of  acetone  is  21.3  which  is  very  large  in  comparison  with  2.2  of  dioxane,  yet 
the  separations  of  the  bands  in  these  are  of  the  same  order  which  denies  the  said 
proportionality  and  hence,  no  such  rotation. 

On  the  other  hand,  the  asymmetrical  double  minimum  offers  a  very  good  explan¬ 
ation  for  the  appearance  of  the  band  pair.  When  a  double  minimum  in  X-H-Cl 
hydrogen-bonded  system  is  assumed,  the  concept  that  the  fundamental  absorption  band 
of  HCl  has  two  components,  is  that  the  unperturbed  ground  level  (i;=0)  of  the  second 
minimum  is  at  a  position  where  it  can  interact  with  the  unperturbed  t;=l  level  of 
the  first  minimum.  The  fact  that  the  intensity  of  the  absorption  band  on  the  long 
wavelength  side  is  weaker  makes  us  infer  that  the  t;'=0  level  of  the  'second  minimum 
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is  lower  than  the  v=l  level  of  the  first  minimum.  If  the  second  minimum  lies  in 
the  region  of  v=2  level  of  the  first  minimum,  the  HCl  fundamental  absorption  can 
not  be  in  two  bands  and  if  the  second  minimum  becomes  still  lower  but  not  equal  to 
the  first  minimum,  the  asymmetry  will  claim  four  absorption  bands  which  are  however 
negated  by  experiment.  The  case  of  the  two  minima  being  equal  in  depth  as  we  see 
on  ammonia  may  be  excepted  from  consideration. 


Fig.  4  illustrates  schematically  a  double  minimum  potential ;  the  level  is  on  arbitary 
scale.  As  a  result  of  perturbation  one  obtains  for  the  perturbed  energy  level  Ea 
and  Eft’*’ 


where  Ei  and  Eo  are  unperturbed  energy  level  values  and  o=Ei—Eo  is  the  separation 
of  the  unperturbed  levels,  tv  is  the  exchange  integral  Hipodr .  The  eigenfunc¬ 
tions  of  the  two  resulting  levels  E^  and  Eb  are  mixtures  of  zero  order  eingenfunctions 
(fi  and  ipo'  corresponding  to  Ei  and  Eo'  respectively.  They  are 


18)  G.  Herzberg,  Spectra  of  Diatomic  Molecules,  P283  (D.  Van.  Nostrand,  New  York). 
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where 


V'4rM;|*+3*+3 

b-\ 

V4|w|*+5»+3 

2V4jM;|>-fo* 

D-\ 

1 

2  V4!m;|*+^ 

and  a*+6*=l . 


(3) 

(4) 


The  vihrationless  eigenfunction  of  HCl  can  be  given  by  where  tpo  is  the  zero 

order  eigenfunction  corresponding  to  the  vibrationless  level  of  the  first  minimum.  The 
matrix  elements  of  transition  moments  of  the  two  bands  v,  and  Vb  are  calculated^’ 

approximately  by  applying  the  conditions  a>b  and  ^tp^* Mipodr  as 

Ma^a{*pi* Mipadv ,  Mb=b[<pi* Mtpock .  (5) 


Therefore  the  intensity  ratio  is  given  by  From  Eqs.  (1)  and  (3),  the  separation 

<5  of  the  unperturbed  levels  and  the  perturbation  energy  |m;|  are  expressed  by 


3={v„-v»)  (I- W)/(l+W) 


{ya—Vb)*—S*  ^ 


(6) 

(7) 


The  relations  between  the  basicity  constant**’  of  proton  accepters  and  values  of  3  and 
|m;|  calculated  by  using  the  values  of  b*/a*  and  ^v=va—vb  are  illustrated  in  Fig.  5. 
Along  with  the  decrease  t)f  proton  attractive  force  of  solvent,  the  values  of  3  show 


basicity 


Fig.  5.  The  relations  between  the  basicity  constant 
and  values  of  3,  \u}\  and  b^ja*. 


19)  W.  Gordy  and  S.  C.  Stanford,  J.  Chem.  Phys.  9  (1941)  204. 
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a  tendency  to  fall  off  as  seen  in  Fig.  5  suggesting  the  rise  of  level  of  the  second 
minimum  relative  to  the  first  and  with  it  the  decrease  in  difference  between  a  and  b 
as  can  be  understood  from  Eq.  (3).  This  means  more  ipx  and  (po  to  mix,  but  the 
hydrogen  bond  becomes  weaker. 

The  experiment  confirms  the  increase  of  intensity  of  the  band  on  the  long  wave¬ 
length  side  due  to  the  intensity  ratio  being  very  much  dependent  on  o  rather 
than  Iwl  as  the  basicity  of  solvent  is  reduced.  |m;1  relates  to  the  height  of  the 
potential  barrier  between  the  two  minima.  Although  this  height  is  diffcult  to  be 
estimated,  it  is  reasonable  to  conclude  that  such  a  barrier  exists  and  such  systems  as 
are  dealt  with  here  do  not  show  a  single  potential  minimum  and  that  the  height  of 
the  second  minimum  is  in  the  vicinity  of  1800  cm“^~1900  cm"*  from  the  vibrationless 
(t;=0)  level  of  the  first  minimum. 

Beside  the  already  discussed  broad  absorption  band,  a  sharp  band  appeared  in  the 
cases  of  diethylether  and  dioxane  at  826  rm~*  and  837  cm"*  respectively  by  the  addition 
of  HCl.  This  band  increases  in  intensity  along  with  HCl  molar  concentration.  It  has 
been  known  that  the  skeletal  vibrational  frequency  shifts  when  molecular  complex  is 
formed.*®’  The  above  mentioned  new  band  is  probably  attributed  to  a  low  frequency 
shift  of  symmetrical  C-O-C  skeletal  frequency  in  the  hydrogen  bond.  Adams  and 
Katz**’  found  on  dioxane  a  similar  band  at  about  840  cm"*  by  the  addition  of  HF. 
This  proves  also  the  formation  of  O-H-Cl  hydrogen  bond.  With  substances  other 
than  diethylether  and  dioxane,  the  frequency  shift  was  not  conspicuous  probably  be¬ 
cause  HCl  molar  concentration  was  insufficient. 

As  to  the  electron  configuration  for  hydrogen  bonding,  the  structure  X*-H:C1" 
should  be  regarded  as  playing  an  important  role.  The  charge  effect  of  non-bonding 
pair  of  electrons  is  connected  with  the  charge  transfer  force  of  Mulliken.  If  the  new 
absorption  band  in  ultraviolet  region  is  due  to  the  charge  transfer  force,  we  should 
consider  it  to  be  the  result  of  transition  between  the  states  N  and  E  corresponding  to 
the  following  two  eigenfunctions 

Vy=CxV[X,  (i/C/)}  -bC,r  {A’+-(//C/)"}  (8a) 

= -C,r  {X,  iHCt)]  -f  CxV  {X*-{HCD-\  (8b) 

where  Ci  and  C,  are  constants  with  the  relation  Ci>Cs.  For  example,  the  highest 
electron-occupied  state  in  the  case  of  acetone  is  the  non-bonding  electron  state  and  its 
ionization  potential  is  9.96  eV,**’  and  the  first  excited  state  of  HCl  molecule  appears 
as  a  continuous  band  starting  at  44()(X)cm"*  and  its  ionization  potential  is  12.7  eV.**’ 

20)  A.  D.  E.  Pullin  and  J.  Me.  Pollock,  Trans.  Faraday  Soc.  54  (1958)  11. 

21)  R.  M.  Adams  and  J.  J.  Katz,  J.  Mol.  Spectr.  1  (1957)  309. 

22)  K.  Watanabe,  J.  Chem,  Phys.  26  (1957)  542. 
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The  exact  position  of  the  maximum  of  this  continuous  band  has  not  been  make  known, 
but  by  analogy  with  the  case  of  HBr  band,  we  may  be  allowed  to  assign  it  above 
50000  cm“^  The  values  of  ionization  potential  given  above  are  for  acetone  and  HCl 
both  in  free  gas  phase ;  in  solutions  they  may  be  different  but  their  relative  positions 
are  as  illustrated  in  Fig.  6.  The  «— ►x*  absorption  band  resulting  from  C=0  of  acetone 
(2790  A  in  hexane  solution)**'**’  shifts  because  of  the  hydrogen  bonding.  But  this  band 
is  different  from  the  7c—*7c*  band**’”'  and  in  most  case  does  not  show  the  “  red  shift  ” 
when  the  hydrogen  bond  is  formed.  Therefore  the  appearence  of  the  band  at  about 
360  m^  owes  likely  to  the  charge  transfer  in  the  molecular  complex. 

XlO^Cm"' 


However,  it  is  a  question  whether  this  can  explain  adequately  the  case  of  absorp¬ 
tion  in  CCU-HCl  solution  superposed  by  the  proper  absorption  of  CCI4.  CCI4  being 
an  inactive  solvent  can  hardly  form  a  molecular  complex  with  HCl  molecule  and  the 

23)  H.  L.  McMurry  and  R.  S.  Mulliken,  Pro.  Nat.  Acad.  Sci.  26  (1940)  312. 

24)  H.L.  McMurry.  J.  Chem.  Phys.  9  (1941)  231,  241. 

25)  S.  E.  Sheppard,  Rev.  Mod.  Phys.  14  (1942)  303. 

26)  H.  McConnel,  J.  Chem.  Phys.  20  (1952)  700. 
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band  in  question  may  in  fact  be  the  band  proper  to  CCI4  or  HCl  but  shifted.  Further, 
it  is  strange  that  the  absorption  bands  of  diethylether  and  isopropylether-HCl  solution 
have  no  maximum.  For  the  present,  as  far  as  the  absorption  band  that  appears  near 
360m/ix  in  the  cases  of  acetone-HCl  etc.  is  concerned,  the  charge  transfer  may  be 
responsisible  for  the  absorption  provided  that  the  band  is  not  caused  by  impurities. 
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